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whi  eh  is  more  than  two  orders  of  magnitude  better  that  the  wavelength 
of  the  microwaves  used. 

The  theory'  of  the  method  has  been  developed  and  the  dominant  features 
are  now  well  understood.  We  have  constructed  an  instrument,  working 
at  10  GHz,  which  has  performed  essentially  in  accordance  with  predictions. 
Using  tills  instrument  we  have  been  able  to  detect  microcracks  in  fatigue 
metal  samples,  with  a width  of  the  order  of  one  micron.  We  were  able 
to  locate  such  cracks  with  an  accuracy'  of  better  than  100  microns.  The 
instrument  has  also  proved  capable  of  detecting  very'  small  changes  in  the 
electrical  properties  of  the  surface  of  a sample.  In  particular,  it  is 
able  to  detect  changes  in  dieleccric  constant  which  arc  substantially  less 
than  1$.  As  another  example,  it  can  measure  the  ferrimagnetic  line  width 
of  a ferrite  sample  - over  an  area  with  dimensions  of  the  order  of  100 
microns.  The  instrument  is  already  capable  of  performing  testing  functions 
which  cannot  be  duplicated  by  other  means.  In  our  view,  it  could  be  further 
developed  for  routine  testing  in  a wide  class  of  applications. 

The  accoustic  surface  wave  technique  differs  from  those  previously  employed 
in  that  we  have  adopted  a holographic  technique,  whereby  it  is  possible 
to  record  the  accoustic  field  along  a single  line,  and  from  this  recover 
the  complete  field.  It.  is  therefore  able  to  detect  defects  in  the  surface 
of  a sample  by  reconst ruction  based  on  a single  record.  Our  method  is 
based  on  the  use  of  a phase  sensitive  laser  probe,  followed  by  subsequent 
computer  reonst ruction.  Using  this  technique  we  have  been  able  to  detect 
very  small  defects  on  a quartz  surface.  It  is  possible  that  the  technique 
can  be  used  effectively  for  determining  surface  flaws  in  highly  polished 
surfaces. 
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ABSTRACT 


We  have  explored  the  feasibility  of  using  two  new  techniques  for  non- 
destructive testing  applications,  one  based  on  the  use  of  microwaves,  the 
other  on  acoustic  surface  waves. 

The  microwaves  scanning  microscope  technique  differs  from  other  microwaves 

i 

n.d.t.  systems  in  that  it  aims  not  only  to  determine  the  existence  of  a 
defect,  but  also  to  locate  it.  The  location  is  effected  with  an  accuracy 
which  is  more  than  two  orders  of  magnitude  better  than  the  wavelength  of 
the  microwaves  used. 

The  theory  of  the  method  has  been  developed  and  the  dominant  features  are 

now  well  understood.  We  have  constructed  an  instrument,  working  at  10  GHz,  i 

which  has  performed  essentially  in  accordance  with  predictions.  Using  this  I 

instrument  we  have  been  able  to  detect  microcracks  in  fatigue  metal  samples, 

with  a width  of  the  order  of  one  micron.  We  were  able  to  locate  such  cracks 

with  an  accuracy  of  better  than  100  microns.  The  instrument  has  also  proved 

capable  ‘of  detecting  very  small  changes  in  the  electrical  properties  of  the 

surface^  a sample.  In  particular,  it  is  able  to  detect  changes  in  dielectric 

constanl/Which .are  substantially  less  than  1%.  As  another  example,  it  can 

measure  the  ferrimagnetic  line  width  of  a ferrite  sample  - over  an  area  with 

dimensions  of  the  order  of  100  microns.  The  instrument  is  already  capable  of 

performing  testing  functions  which  cannot  be  duplicated  by  other  means.  In 

our  view,  it  could  be  further  developed  for  routine  testing  in  a wide  class 

of  applications. 

The  acoustic  surface  wave  technique  differs  from  those  previously  employed 
in  that  we  have  adopted  a holographic  technique,  whereby  it  is  possible  to 
record  the  acoustic  field  along  a single  line,  and  from  this  recover  the  complete 
field.  It  is  therefore  able  to  detect  defects  in  the  surface  of  a sample  by 
reconstruction  based  on  a single  record.  Our  method  is  based  on  the  use  of 
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a phase  sensitive  laser  probe,  followed  by  subsequent  computer  reconstruction. 
Using  this  technique  we  have  been  able  to  detect  very  small  defects  on  a 
quartz  surface.  It  is  possible  that  the  technique  can  be  used  effectively 
for  determining  surface  flaws  in  highly  polished  surfaces. 
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SECTION  A.  MICROWAVE  SCANNING  TECHNIQUE 
1.  INTRODUCTION 

This  is  the  Final  Report  covering  the  progress  in  the  last  quarter  from 
1 October  to  31  December,  as  well  as  the  work  done  throughout  the  project.  The 
work  on  the  microwave  scanning  technique  is  discussed  in  section  A,  whilst  that 
on  acoustic  surface  wave  holography  is  presented  in  section  B. 

We  believe  that  both  of  the  techniques  have  been  shown  capable  of  performing 
certain  non-destructive  testing  functions  which  go  beyond  the  capabilities  of 
existing  methods.  A brief  overall  assessment  of  the  program  and  plans  for  future 
work  are  presented  in  section  C. 

Experiments  have  been  carried  out  to  assess  the 
potential  use  of  the  microscope  in  non-destructive  testing 
applications.  These  demonstrate  that  the  microscope  is 
capable  of  resolving  object  details  separated  by  distances 
more  than  tuo  orders  of  magnitude  smaller  than  the  illuminat- 
ing wavelength.  The  sensitivity  of  the  instrument  to  changes 
in  surface  topography  of  metals  has  been  exploited  to  detect 
micron  wide  fatigue  cracks  in  metals;  the  microscope  is  also 
very  sensitive  for  the  detection  of  breaks  in  thin  metallic 
films  on  substrates,  this  makes  it  useful  in  the  testing  for 
cracks  in  thin  film  circuitry. 

A technique  whereby  we  can  selectively  image  small 
changes  in  electric  and  magnetic  properties  of  the  object  is 
described.  The  limited  depth  of  field  of  the  microscope 
enables  us  to  characterise  the  properties  of  both  dielectric 
and  ferrite  thin  films.  It  would  be  extraordinarily  difficult 
to  visualise  these  changes,  non-destructi vely , in  any  other 


way . 
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2.  BASIC  PRINCIPLE  OF  METHOD 

The  basic  principle  of  the  method  depends  on  the  use  of  an  aperture  in  a 
metallic  plate  which  is  small  compared  with  the  illuminating  wavelength*  Fig. 1.1. 
We  view  the  object  through  the  aperture.  If  the  object  is  transparent,  one  can 
detect  the  transmitted  wave  through  the  hole;  if  it  is  reflecting,  one  can 
detect  the  scattered  wave  through  the  aperture.  By  scanning  the  sample  under 
the  pupil  and  recording  the  signal  we  can  eventually  scan  the  whole  sample  with 
a resolution  which  is  somewhat  less  than  that  of  the  aperture  dimension  but  more 
than  two  order  of  magnitude  less  than  the  wavelength. 


The  "objective"  of  our  microscope  consists  of  an 
aperature  in  one  of  the  side  walls  of  a resonant  cavity, 
which  in  our  case  is  a microstrip  resonator  (Fig,  1,2;, 

The  object  beneath  the  hole  perturbs  the  frequency  of  the 
resonator.  By  operating  on  the  linear  portion  of  the 
resonator  Q curve,  we  can  translate  this  frequency  shift  v. 

into  a corresponding  amplitude  change,  which  is  demodulated 
in  the  receiver  and  displayed  on  a X-Y  plotter.  The  object 


is  illuminated  by  th9  evanescent  fields  associated  with  the 
small  aperature,  Ue  must  therefore  detect  a signal  which  is 
very  small  compared  with  the  illumination  power  in  order  to 
obtain  a significant  improvement  over  the  classical  resolution. 
Moreover,  this  signal  is  superimposed  on  a much  larger  back- 
ground radiation  so  that  the  contrast  will  be  poor.  This 
deficiency  can  be  overcome  by  vibrating  the  object  at  an 
acoustic  frequency.  The  required  signal  is  then  "tagged" 
by  this  modulation  frequency,  and  is  readily  separated  in 
the  receiver  system. 


Illuminating  Detector  for 
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In  the  case  of  some  materials  it  is  possible  to  apply  the  modulation  in  an 
electrical  rather  than  a mechanical  form.  For  example,  in  the  case  of  a ferrite 
object  it  is  possible  to  modulate  the  magnetic  field,  by  modulating  a current 
in  a suitable  auxiliary  coil. 

3.  THEORY  OF  MICROSCOPE 

In  aperture  microscopy,  just  as  in  ordinary  microscopy,  we  are  interested 
in  assessing  the  resolution  and  sensitivity  of  the  instrument  to  various  materials 
and  material  defects.  In  this  section,  we  shall  present  approximate  theories  to 
explain  the  .-esponse  of  the  instrument  to  both  simple  variations  in  the  surface 

topography  as  well  as  changes  in  the  electric  and  magnetic  permittivity  of  the 
sample. 


Tn  order  to  predict  the  performance  of  the  system, 
we  have  used  an  approximate  theory  to  assess  the  signal 
strength  as  a function  of  position  of  a metallic  plane 
belou  the  aperture.  However,  in  practice,  ue  are  far 
more  interested  in  how  small  deviations  in  the  flat  plane 
affect  the  signal,  Ue  have  therefore  developed  a first 
order  theory  to  assess  the  perturbation  in  signal  from  a 
defect,  simulated  as  a slot  in  a metal  object. 

In  the  use  of  apertures  which  are  small  as  compared 
uith  one  wavelength,  it  is  meaningful  to  distinguish  between 
electric  and  magnetic  contrast;  this  is  used  to  explain  the 
ability  of  the  microscope  to  image  smail  perturbations  in  the 
relative  dielectric  constant  as  well  as  variations  in  the 
permeability  of  the  object. 


4. 


3.1  Aperture  Superresolution 

The  imaging  properties  of  the  microscope  are  determined  by  the  evanescent 
fields  produced  by  the  cavity  aperture.  Our  ability  to  obtain  high  resolution 
is  based  on  the  fact  tuat  this  evanescent  field  contains  spatial  frequencies 
which  are  much  higher  than  the  inverse  wavelength.  We  have  looked  at  this 
issue  quantitatively  by,  computing  the  Fourier  transforms  of  approximate 
versions  of  the  aperture  field  distribution.  v 

Using  the  magnetic  dipole  approximation  to  the 
aperture  field,  ue  computed  the  transverse  spectral  density 
function  in  the  X direction^ ( due  to  the  symmetry  of  the  fields 
this  is  the  same  as  the  spectral  density  in  the  Y direction), 
the  result  is  shcun  in  Fig.  2.5a,  The  calculation  was 
performed  for  an  aperture  diameter  ue  have  been  using  in 
experiments  of  "'VlQO,  and  for  a typical  operating  distance 
of  Z *=  ^/1000  auay  from  the  diaphragm.  The  shaded  region  v. 

in  the  figure  shows  that  the  spatial  frequencies  which  are 
mors  than  two  orders  of  magnitude  higher  than  the  inverse 
wavelength  are  transferred  by  the  aperture.  It  is  precisely 
this  part  of  the  spectrum  on  which  we  rely  to  obtain  super- 
resolution  imaging  of  the  object. 

It  was  of  interest  to  evaluate  the  spectral  power 
density  in  the  Z direction.  This  was  computed  from  Bethe's 
approximation  to  the  component  and  is  shown  in  Fig.  2,5b.  ^ 

The  figure  indicates  that  we  in  fact  have  super-resolution 
in  depth  as  well  as  in  the  transverse  directiun.  It  also 


shows  that  the  resolution  decreases  very  rapidly  as  we  move 
away  from  the  aperture  in  the  Z direction.  This  result  has 
tuo  important  implications.  First,  as  any  practical  instrument 
will  have  a limited  dynamic  range,  we  must  operate  as  close 
to  the  aperture  as  possible  if  ue  are  to  detect  these  high 
spatial  frequencies  to  obtain  the  high  resolutions  predicted. 
Secondly,  the  fact  that  the  depth  of  field  of  the  microscope 
is  itself  less  than  an  aperture  diameter,  suggests  that  ue 
might  see  surface  changes  which  would  be  extraordinarily 
difficult  to  detect  in  any  other  way.  In  Section  7.3  we 
shall  present  the  results  obtained  on  the  imaging  of  varia- 


tions in  the  properties  of  thin  films  which  confirms  this 
important  property. 

The  spectral  density  functions  considered  are  useful 
to  illustrate  the  principle  of  aperture  supBr-resolution ; the 
results  indicate  that  the  aperture  facilitates  the  transfer 
of  very  high  spatial  frequencies,  so  that  one  might  expect 
to  achieve  resolutions  much  less  than  an  aperture  diameter. 

In  practice,  the  resolution  ue  can  achieve  depends  on  the 
interaction  of  the  aperture  fields  with  the  particular  object 
under  investigation  as  well  as  the  overall  dynamic  range  of 
the  microscope.  The  determination  of  the  object/aperture 
field  interaction  requires  a detailed  field  calculation  for 
each  particular  type  of  object;  consequently,  we  have  not 
attempted  to  solve  these  very  complex  problems.  However, 
we  can  say  that  for  most  objects,  the  attainable  resolution 
will  be  of  the  order  of  the  aperture  diameter  ( tylOO  in  our 

case) . 


1 
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The  basic  concept  of  the  microscops  is  shown  in 
fig.  1.2.  Ue  illuminate  a thin  metallic  screen  with  an 
aperature  of  diameter,  a,  which  is  very  much  sr.,aller  than 
the  illuminating  wavelength,  X,  , This  small  aperature  gives 
rise  to  a distribution  of  electromagnetic  waves  which  are 
evanescent  in  a direction  normal  to  the  diaphragm.  Ue 
must  therefore  detect  a signal  which  is  very  small  as  compared 
to  the  illuminating  power.  Moreover,  this  signal  i3  super- 
imposed on  a much  larger  background  radiation,  so  that  the 
•contrast'  will  be  very  poor.  The  loss  of  signal  may  always 
be  recovered  by  the  use  of  sufficiently  long  integration  times. 
The  more  serious  problem  of  the  loss  of  contrast  can  be  over- 
come by  vibrating  the  object  at  an  audio  frequency,  and 
subsequently  searching  for  this  modulated  component  in  the 
receiver « 


Bethe  has  shown,  that  fields  on  the  object  side  of  the  diaphragm  (Z>0)  can  be 
approximated  by  a pair  of  magnetic  and  electric  dipoles.  If  the  screen  is 
illuminated  so  as  to  produce  a predominantly  magnetic  dipole 
orientated  in  the  plane  of  the  screen,  then  the  ob ject/aperatur 
field  interaction  gives  rise  to  a dipole  field  in  the  region 
Z<-0  a^ove  the  diaphragm.  The  strength  and  orientation  of 
this  field  will  contain  the  desired  information  op  the  object, 
for  the  case  when  the  electric  dipole  predominates  and  is 
orientated  in  a plane  normal  to  the  aperture,  only  the  strength 
of  the  dipole  fields  in  a region  Z < 0 will  provide  information 
on  the  properties  of  the  object. 
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An  alternative  way  of  looking  at  the  theory  of 
operation  is  to  consider  a technique  which  applies  more 
closely  tc  the  theory  which  will  be  developed  in  the 
following  sections.  In  our  case,  the  dipole  radiates  into 

f 

j 

j a microstrip  resonator.  The  object  under  the  aperture  will 

I 

perturb  the  resonant  frequency  of  the  resonator  as  well 

V. 

changing  the  unperturbed  Q of  the  resonator.  It  uill  be 

J 

shown  in  Section  3.2  that  we  are  primarily  interested  in 
changes  in  the  resonant  frequency.  By  operating  on  the 
linear  portion  of  the  resonator  Q curve,  we  can  translate 
this  frequency  3hift  to  a corresponding  amplitude  change, 
unich  is  then  demodulated  in  the  receiver  system  and  displayed 

The  sources  of  object  contrast  ere  found  in  the 
distribution  of  electric  and  magnetic  permittivities  and  in 
the  distribution  of  conductivity  which  alter  the  energy  stored 
in  the  vicinity  of  the  aperture.  If  there  are  sufficient 
gradients  of  these  quantities  within  the  dipole  field,  there 
uill  also  be  a contribution  to  the  signal  arising  from  the 
distortion  of  thb  dipole  field. 

3 . 3 Effect  of  Metallic  Plane 

I 

! In  order  to  obtain  some  estimate  of  the  signal 

i strength,  we  have  calculated  tne  frequency  perturbation 

f 

j caused  by  e metal  object  which  is  vibrated  so  as  to  close 

the  aperture  at  one  part  of  the  cycle,  and  which  is  removed, 

| effectively  outside  the  range  of  the  dipole  near- field, 

I half  a cycle  later. 
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Using  simple  perturbation  theory,  assuming  that  the  change  in  resonant 
frequency  Sf  is  given  by 

ME , 

2 

whereU^  & are  the  electric  and  magnetic  energies  stored  in  the  near  field  of 
the  aperture,  and  is  the  total  energy  stored  in  the  resonator.  For  the  case 
of  an  open  resonator^,  the  signal  power,  modulated  at  the  frequency  fm  was 
calculated  at  the  optimum  frequency  (displaced  from  the  resonent  frequency  by 
a fraction  3^/Q,  took  the 

S'  0 

_2 

kQ  = 2jt/\  is  the  wave  number.  Typically***  10  and  the  aperture  diameter 
a = A/100,  wo  being  the  minimum  spot  diameter;  the  signal  power  is  approximately 
52dB  below  the  input  power. 

Ue  have  evaluated  the  signal  amplitude  uhich  ue 
obtain  uhen  a metallic  object  is  removed  from  contact  with 
the  aperture  to  a position  outside  the  range  of  the  dipole 
field.  However,  ue  need  finer  grained  information  - how 
the  signal  varies  with  the  amplitude  of  vibration  and  with 
mean  position  of  the  object,  Essentially  ue  are  concerned 
with  a quasi-static  field  theory.  If  ue  can  regard  the 
object  as  being  a good  conductor,  its  effect  can  be  taken 
into  account  using  image  theory.  The  image  of  the  aperture 

is  itself  imaged  in  the  aporture  plane  and  this  too  can  be 
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incorporated  in  the  theory  (Tig.  3.3),  the  process  can  be 
continued  in  the  'hall  of  mirrors”,  the  effect  of  each 
successive  i^age  to  the  total  field  being  increasingly 
less  important.  As  long  as  the  distance  between  the  object 
and  the  image  plane  is  large  as  compared  with  the  aperture 
diameter,  we  have  a formulation  which  is  very  nearly  rigorous. 
In  practice  we  must  daal  with  object/aperture  separations 
which  are  comparable  to  the  aperture  diameter.  In  this 
regime,  our  theory  becomes  progressively  less  reliable. 
Nevertheless,  it  should  grasp  the  most  important  features 
of  the  Situation  and  provide  a useful  guide  to  design. 

The  calculation  of  the  signal  strength  involves 
evaluating  the  seif  energy  of  the  dipole  and  its  FI  images, 
as  well  as  the  cross  or  mutual  energy  terms  associated  with 
the  dipole  and  each  of  its  FI  images. 

The  basic  t/chniqua  we  have  adopted  begins  by 
approximating  the  aperture  fields.  The  object  is  assumed 
to  be  a perfect  conductor  and  its  effect  may  then  be  taken 
into  account  using  image  theory.  The  resultant  fields  in 
the  vicinity  of  the  aperture  are  then  represented  as  a set 
of  plane  waves.  The  change  in  energy  stored  as  the  screen 
is  displaced  is  then  computed  in  the  frequency  domain  rather 
than  in  the  configuration  domain.  In  the  first  part  of  the 
analysis  we  deal  with  the  one  dimensional  problem  - the  case 
of  the  infinite  slot  aperature  (Fig  3.4);  this  is  then 
extended  to  include  the  two  dimensional  problem  of  a square 
aperture . 

First,  ue  assume  that  the  function  g (x,0)  is  a 
good  approximation  of  the  aperture  field  distribution  in  the 
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Figure  3.3  Imaging  of  Dipole  in  Metallic 
Object  and  in  Conducting  Screen 
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The  field  in  the  vicinity  of  tne  hole  may  be 


represented  by  a set  of  plane  uaves: 

fOO  ♦ I 

3(x'0)= 


..(3.3  ) 


- o O 


where  G(fx)  is  the  fourier  transform  of  g(x,0).  Each 
term  of  equation  3.18  is  a plane  U3ve,  and  ue  may  express 
the  field  distribution  at  a position  2 away  from  the  aperture 

as: 


9(K'*)= 


r°  -jkx  -jk.* 
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(f*)e 


A(f») 


-oO 


..(3. 4 ) 


From  the  uave  aquation, 


kz  is  given  by: 


..(3. 5 ) 


Ue  are  operating  in  the  very  near  field  of  the  aperture, 
i.e.  in  the  induction  zone,  and  consequently  for  the  most 

o o 

important  field  components  the  approximation  k >>k 

X 0 

will  be  valid ; and  3,20  may  be  rewritten  as: 


..(3.  6 ) 


From  3, A and  3,6  ue  obtain: 


g(x,i)  t 
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The  energy  stored  per  unit  length  as  the  screen  is  displaced 
from  2 = 0 to  2 = d,  for  the  source  c?.i  be  shown  to  ber 
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A method  using  Discrete  Fourier  Transforms  is  particularly 
useful  in  evaluating  this  equation  os  the  technique  readiiy 
lends  it'elf  to  computer  numerical  methods. 

From  equation  3.  8,  the  energy  per  unit  length  of 
the  slot  is: 

/.  x.  , , ra  -ik 

£ | G ( f n ) | (A{) 


V/  = 60L 


n-0 


<U 


..(3.  9 ) 


As  kn  » 2TTfn  oncJ  = nA  ue  maY  rewrite  equation  3.®'  as: 


v / « £oL  L 1 

4TT  n = o J 


..(3.10) 


Let  us  assume  that  the  aperture  field  distribution 
g(x,0),  is  a gaussian  given  by: 


3 (x.  o)  =■  e 
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..(3.  11) 


^L!W!W4^'--.S!A','J  -^rr^ 


15 


m 


I 


15. 


Xq  being  the  gaussian  decay  constant* 

If  the  effect  of  the  conducting  planeis  taken  into  account 
using  image  theory  and  using  equation  3.10  and  3.11  a typical 
component  for  the  M'th  image  is: 


«(*,*)  =lbtl  G(>) 

J x ft -0 


A rfr  \ -jKnX-kn(^f^m'0 


..(3. *2) 


•where : ‘ d(o)  refers  to  the  case  of  the  source  alone 

d(m)  » md;  if  m is  even 
and  d(m)  * (m  + 1)  ; if  m is  odd 

The  self  energy  terms  for  the  source  and  its  N images  becomes: 


I 


. . i//fe^V!ffl4rln 

VLf=M?  7 7 1^7 — le 


-qrr 


Z,  L 

w*o  n-  o 


l 


- e 


£(AM+<1)~ 


,.(3.  13) 


is 


I 

8 

£ 

£ 

if 


i 

t 

t 

f? 

i 

! 

!<' 

t 

i 

P 


16. 


while  the  cross  or  mutual  energy  terms  are  given  by: 


M M W 
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, = £5X0 

^OSS  A Li  LtU 

tn-o  S-nH  n-o 


-ap^A(rr.)U(s)j 

— ' S j 


..(3.14  5 


where,  in  this  case:  d(m)  = (1  + m/2)d  if  m is  even, 

and  d(m)  = (1  + m)d/2  if  m is  odd 

These  results  may  be  extended  to  the  tuo  dimensional 
case  of *the  square  aperture  if  the  following  assumptions  are 
made:  (a)  we  assume  that  the  aperture  field  distribution 

g(x,y,0)  is  separable,  this  allows  the  complicated  tuo 
dimensional  manipulations  to  be  reduced  to  the  more  simple 
one  dimensional  manipulations;  (b)  we  further  assume  that 
the  separable  functions  are  both  symmetric.  uie  may  then 

write  equation  3.12  for  the  tuo  dimensional  case  as: 


H Nf  -3(knX*kfU)-  (fy\ 


• •(3,15  ) 


In  the  above,  we  have  assumed  that  g(x)  is  a gaussian  as 
before  and  g(x)  = g(y),  ue  may  therefore  put  G ( f n ) = G(fp); 
the  subscript  ''p"  referring  to  the  summation  performed  in 
the  Y plane. 

It  may  be  shown  that  the  self  enCigy  for  the  square 


aperture  is: 
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and  the  mutual  energy  terms: 


MM  N N 
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The  total  energy  stored  in  the  quasi-static  aperture  fir'.d 
is  simply: 


total 


U + U , „ 

cross  self 


The  relative  decay  of  the  signal  as  the  screen  is  displaced 
is  simply  the  differential  of  the  total  energy  stored  with 
respect  to  the  movement,  d,  from  the  aperture  in  the  Z 
direction : 


total 


z = d 


A computer  programme  uas  developed  to  calculate  the  signal 
as  a function  of  screen  movement  for  the  source  and  its 
multiple  images.  The  results  are  shown  in  fig.  3.5,  where 
it  must  be  remembered  that  by  the  time  the  object  has 
approached  to  within  one  aperture  diameter,  the  theory  becomes 
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3.4  Perturbation  Due  to  Slot  in  Metal  Object 

In  practice,  we  are  concerned  with  how  small  deviations  in  the  surface 
topography  and  material  changes  affect  the  signal.  Even  for  the  simplest 

defect  geometry,  such  as  a slot  in  a metal  object,  an  electro- 
static field  solution  would  be  very  complex  indeed.  Ue  have 
therefore  made  no  attempt  to  solve  the  problem  rigorously. 
Instead,  we  have  developed  an  approximate  theory  uhich  provides 
a first  order  guide  in  assessing  the  perturbation  in  signal 
arising  from  a defect,  simulated  as  a slot  in  the  metal  object 
surface . 

Ue  wish  to  evaluate  the  change  in  energy  stored  in 
the  vicinity  of  the  aperture  by  the  presence  of  a slot  of 
of  width,  U,  depth,  D,  and  length,  L,  located  symmetrically 
below  a square  aperture  in  a metallic  object  (Fig.  3.6), 

The  slot  is  represented  as  a waveguide  and  the  field  components 
in  the  slot  are  calculated.  The  energy  stored  in  the  slot 
can  then  be  determined.  There  will  also  be  a contribution 
to  the  energy  stored  by  the  slot  as  a result  of  the  perturba- 
tion of  the  fields  above  the  slot.  This  contribution  has 
been  neglected  as  we  only  require  an  order  of  magnitude 
estimate  of  the  perturbation  due  to  the  slot,  Ue  have  also 

neglected  second  order  effects,  suco  as  the  interaction  of 
the  eJye  of  the  slot  with  the  side  of  the  aperture. 

Let  us  consider  the  case  when  tho  aperture  fields 
may  be  regarded  as  due  to  those  of  magnetic  dipole.  This 
corresponds  to  exciting  the  slot  uith  a magnetic  field  com- 
ponent H^,  the  electric  field  component  being  zero. 
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The  field  components  in  the  slot  were  evaluated  and  used  to 
calculate  the  energy  stored  in  the  slot: 


Equation  3,£g  indicates  that  for  thin  deep  slots  (to«J)), 
the  signal  related  to  the  slot  depth  uill  saturate. 
Experimental  results  which  verify  the  thenr*,  uill  be 
presented  in  Chapter  6.  The  equation  also  predicts  that 
the  signal  should  be  proportional  to  the  square  of  the  uidth 
of  the  slot.  Experimental  results  indicate  that  the 
dependence  of  signal  on  uidth  is  in  fact  more  nearly  linear. 

The  theory  accurately  predicts  the  decay  of  signal 
with  the  depth  of  the  slot,  but  not  the  variation  of  signal 
with  the  uidth  of  the  s’nt.  This  apparent  anomaly  can  be 
explained  if  ue  consider  hou  each  quantity  arises  in  the 

final  expression.  The  decay  of  the  signal  is  related 


primarily  to  the  attenuation  constant  in  the  slot  and  this 
parameter  can  be  calculated  accurately.  However,  the  change 
in  signal  with  uidth  of  slot  is  determined  by  the  precise 
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field  configuration  inside  and  outside  the  slot.  Us  have 
not  determined  these  fields  with  any  great  accuracy,  thus 
the  discrepancy  in  the  theoretical  and  experimental  results. 

For  slots  which  have  U«D,  the  energy  stored  in  the 
slot  may  be  approximated  to: 


1) 

UftlOT 


IL  K\tJ\ 
/* 


..(3  21  ) 


In  order  to  evaluate  the  perturbation  in  signal 
due  to  the  slot  in  a metallic  object  positioned  at  a distance, 
d,  auay  from  the  aperture,  we  must  first  calculate  the  energy 
stored  due  to  the  aperture  and  plane  in  the  absence  of  the 
slot.  The  energy  stored  in  the  near  field  of  a square 
aperture  in  the  presence  of  a metallic  plane  was  calculated 
in  Chapter  3.2.  The  energy  stored  'or  the  dipole  and  its 
first  image  uas: 

Let  us  consider  a slot  uhich  has  U D ~ this 
corresponds  to  the  type  of  fatigue  cracks  ue  have  been 
detecting.  From  equations  3,44  and  3.45  the  perturbation 
in  energy,  A U due  to  the  slot  in  a metallic  object  positioned, 
d,  auay  from  ths  aperture  plane  is: 


0'393 


..(3.22  ) 
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If  a unity  signal  to  noise  ratio  is  used  as  the 

criterion,  the  minimum  detectable  perturbation  of  the  system 

is  Au  where: 
min* 


A0mm  - ytVsl)  ..(3.24) 

For  a typical  height  of  d = a/10  of  the  metallic  plane  below 
the  aperture,  of  diameter  a = 300  yt  and  length  of  slot 
L »=  300  y,  the  predicted  minimum  slot  width  which  can  be 
detected  for  a system  s/N  ratio  of  45  dB  is  about  4 y. 
Experimental  results  (Chapter  6)  indicate  that  equation  3.47 
does  predict  the  order  of  magnitude  of  the  detection  sensit- 
ivity. 
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3.5  Dielectric  and  Ferrite  Imaging 

We  have  mentioned  that  the  microscope  is  sensitive 
to  changes  in  the  distributj c.i  of  the  electric  and  magnetic 
permittivity  of  the  object.  In  this  section,  we  consider 
a scheme  which  uill  enable  us  to  utilise  the  microscope 
selectively  to  image  variations  in  the  relative  dielectric 
constant  and  variations  in  the  permeability  of  the  object. 

The  imaging  properties  of  the  microscope  may  be 
explained  by  studying  the  changes  in  the  energy  stored  in 
the  vicinity  of  the  aperture  by  the  object.  To  a first 


25. 


approximation,  the  aperture  fields  can  be  represented  by 
either  a magnetic  or  an  electric  dipole  depending  on  the 
form  of  illumination.  Let  us  consider  the  energy  stored 
in  the  induction  field  of  the  aperture  for  both  the  electric 
and  the  magnetic  dipole  case.  It  Can  be  shown 
that  for  the  magnetic  dipole,  the  principle  energy  storage 
mechanism  is  magnetic  and 


£/„.  - O-o^Hoa 


While  for  the  electric  dipole  it  is  the  electric  energy 
stored  which  is  of  consequenca: 


na 


* O-Oll  at  £° 


The  above  equations  only  apply  when  the  material  in  the  half 
space  below  the  aperture  is  filled  with  a material  of 
dielectric  constant,  S. , in  the  electric  dipole  case,  and  a 
material  of  permeability,  j*.  , in  the  magnetic  dipole  case 
If  the  perturbation  in  either,  £.  , or,  p,  is  small, 


(7mh  06 

UelU  °C  £ 


Houever,  if  there  are  significant  gradients  in  the  electric 
or  magnetic  permittivity  of  the  object,  tnere  will  also  be 


a contribution  to  the  signal  due  to  the  distortion  of  the 
aperture  fields,  and  the  signal  will  no  longer  be  proportional 
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to  these  quantities  i.e.  first  order  perturbation  theory 
fails. 

Tuo  important  points  emerge  from  this  very  approx- 
imate theory.  If  ue  suitably  illuminate  the  aperture  to 
obtain  a predominantly  electric  dipole,  the  microscope  uill 
be  very  much  more  sensitive  to  changes  in  the  relative 
dielectric  constant  than  to  changes  in  the  permeability 
of  the  object.  Similarly,  for  an  essentially  magnetic 
dipole  aperture  field,  the  instrument  will  be  more  sensitive 
to  variations  in  the  permeability  of  the  sample.  Experiments 
to  corroborate  these  findings  are  presented  in  Chapter  7,3. 
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4.  RESONATOR  AND  APERTURE  DESIGN 

i 

The  "heart"  of  a microscope  is  the  objective;  i 

in  our  case,  this  consists  of  a resonator  with  an  aperture  j 

j 

in  one  of  the  side  walls.  In  this  chapter  ue  shall  deal 
with  the  design  and  optimisation  of  the  microscope  "objective". 

The  main  design  criterion  is  to  maximise  the 

object/aperture  field  interaction.  It  was  shown  in  Section 

3,  that  the  signal  was  proportional  to  the  energy  stored 

in  the  fields  outside  the  aperture,  Uo,  and  inversely 

proportional  to  the  energy  stored  in  the  resonator,  Ur» 

Ue  therefore  need  to  increase  the  ratio  U /u  The  increase 

o'  r. 

in  UQ  may  be  effected  by  confining  the  fields  illuminating 
the  aperture  into  the  smallest  possible  area  so  as  to  have 
a maximum  field  intensity  at  the  aperture. 


The  choice  of  the  type  of  microwave  resonator  other  than  optimising 
the  signal  to  power  ratio,  depends  on  factors  such  as  the  size,  ease  of  fabrication, 
and  stability  of  the  resonator.  These  considerations  indicated  that  a micro- 
strip resonator  would  have  significant  advantages  ove^  the  open  resonator  as  the 

microstrip  objective. 

4.1  The  Microstr.ip  Resonator 

Several  forms  of  microstrip  resonators  were  invest- 
igated, both  experimentally  and  theoretically  with  the  aim  of 
optimising  the  detection  signal.  They  are  shown  in  fig.  4.1. 

The  investigation  revealed  that  the  linear  resonator  was  able  to  confine  the 
substrate  fields  to  a much  smaller  area  than  the  disc  re$on?tor  ^ ),  or 
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Having  chosen  the  linear  resonator,  we  pro- 
ceeded to  select  the  optimum  design  parameters.  In  the 
following  analysis  ue  shall  consider  a linear  /\s  open  circuit 
resonator  for  convenience;  the  choice  of  a rather  than 
the  optimum  \/2  resonator  is  explained  in  the  next  section. 


The  linear  resonator  we  are  analysing  is  shown  in 
Fig.  4.3.  It  consists  essentially  of  two  discontinuities 
in  the  microstrip  line,  of  length  ^5 . We  shall  assume  that 
the  gaps  in  the  line  have  been  adjusted  for  optimum  coupling, 
and  that  the  aperture  is  positionad  in  the  ground  plane  so 
that  ue  have  a maximum  elactric  field  at  the  aperture.  In 
this  case,  the  aperture  fields  are  essentially  those  produces 
by  an  equivalent  electric  dipcle.  Fig.  4.2  shows  the  field 
configuration  for  the  case  of  a strip  conductor  above  a 
conducting  plane  separated  by  a dielectric  sheet.  We  shall 
also  assume  that  the  propogating  microstrip  ir.c^e  can  be 

(3  ) 

approximated  closely  by  a 7EM  mode'  . 


It  can- be  shown  that  the  signal  to 

input  power  is  related  by: 


0-16  ’6  Q 
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While  the  energy  stored  outside  the  aperture  for  an  slectric 
dipole, 


. * 

Using  this  analysis  we  have  calculated  the  S/P  ratio  to  be:- 


* (The  detailed  analysis  may  be  found  in A.Husain’sphO  Thesis). 
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uhere  £’  is  the  electric  field  normal  to  the  plane  of  the 
aperture.  It  now  remains  to  evaluate  the  energy  stored 
in  the  resonator  and  determine  the  field  incident  on  the 
iperture  in  the  ground  plane. 

A detailed  theory  to  evaluate  E’  and  Ur  uould  be 

very  complicated.  As  ue  only  required  an  approximate  guide 

to  the  optimum  resonator  design  parameters,  ue  have  developed 

(3  ) 

a first  order  theory  based  on  Wheeler’s  analysis'  . 

Thus  for  a particular  aperture  size  and  frequency 
of  oparation,  ue  must  select  the  parameters:  h,  u/h, 

and  0 in  order  to  optimise  the  5/P  ratio.  Equation  4.14 
suggests  that  for  large  S/P  ue  need  to  increase.  Q and  u/h, 
uhile  decreasing  £ and  h.  Houever,  these  four  quantities 
are  not  simply  related.  For  a given  h,  a decrease  in  £.  r 
increases  the  characteristic  impedance  ZQ  and  hence  changes 
the  Q.  The  radiation  that  occurs  from  an  open  circuit 
microstrip  resonator  also  varies  with  £ , ZQ,  h,  and  operat- 

ing  frequency.  Troughton'  ' has  shoun  that  if  these  factors 
are  taken  into  account  and  ailouance  made  for  the  resistive 
losses,  then  the  Q (for  a given  h)  is  a maximum  at  a particular 
impedance. 

In  the  numerical  optimisation  procedure,  the  value 

of  ZQ  for  a particular,  u/h,  uas  obtained  from  the  uork  of 

Yamashita  and  Mittra^  \ and  the  corresponding  Q from  the 

f 5 ) 

theoretical  uork  of  Troughton'  '.  The  effective  dielectric 
constant  uas  calculated  from  equation  4.6.  Using  equation 


4,14,  it  can  be  shoun  that  for  h = 0.25  mm,  the  signal  to 


h * 0.625  mm  or  1.25  mm;  consequently,  h,  was  set  to  the 
minimum  value  (0.25  mm)  for  which  experimental  results  were 
available. 

fig.  4.5  shows  that  the  S/P  varies  with  the  parameter 
w/!i  for  different  values  of  k.  The  S/P  ratio  was  calculated 
for  an  operating  frequency  of  10  GHz,  and  an  aperture  diameter 
of  300  y,  It  is  clear  from  the  figure  that,  subject  to 
assumptions  made  earlier,  the  optimum  parameters  of  the 
linear  resonator  are: 

h * 0.25  mm 

2w/h  = 1.25  ± 0.25 

k « 9-1 

ZQ  3 45  i 4 ohms 

The  variation  in  the  various  resonator  parameters  shown  above 
were  calculated  for  a 5%  reduction  in  the  signal  to  power  ratio 
The  S/P  ratio  for  the  above  parameters  were  calculated  to  be 
0.5  x 10"4. 

This  effect  was  taken  into  consideration  in  the  calculation 
of  the_  optimum  resonator  parameters  using  Wheelers  theory.  The 

result  of  the  correction  is  shown  in  Fig,  4,7  where  we  have 
plotted  the  optimum  S/P  vs.  u/h  curve  (for  k = 8,9),  as  well 
as  the  uncorrected  version  for  comparision.  The  effect  of 
this  correction  is  to  alter  the  optimum  value  of  u/h  slightly 
from  w/h  = 1.25  to  u/h  = 1.4. 
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Having  calculated  the  optimum  design  parameters,  ue 
can  now  discuss  some  of  the  technological  aspects  of  resonator 
construction.  The  microstrip  resonator  uas  fabricated  on  a 
copper  coated  alumina  substrate  using  conventional  photolith- 
ographic techniques.  The  substrates  used  had  a surface  finish 

• In  practice,  ue  found  that  the  optimum  0.25  mm 
thick  substrates  uere  very  fragile,  and  had  a tendency  to 
crack  easily  under  slight  pressure  (such  as  an  accidental 
collision  uith  the  vibrating  object).  Ue  therefore  resorted 
to  use  a thicker  substrate  of  h = 0.5  mm,  even  though  this 
meant  ue  had  to  sacrifice  a loss  in  signal  of  about  3dB. 

In  our  experiments,  the  microstrip  circuit  uas 
mounted  on  a jig.  This  uas  constructed  to  allou  us  freedom 
to  scan  large  specimens  beneath  the  aperture  plane;  uhich 
required  the  probes  to  be  arranged  to  couple  to  the  microstrip 
at  an  angle  (Fig.  4.8);  resulti”n  in  some  radiation  from  the 
bends  in  the  probes.  There  uas  also  radiation  from  the  edge 
of  the  microstrip  circuit.  One  consequence  of  this  uas  that 
a spurious  signal  uas  generated;  this  arose  from  the  modulation 
of  the  energy  leaking  from  the  substrate  by  objects  (vibrating 
at  the  same  acoustic  signal  frequency)  outside  the  resonator. 

Ue  found  that  this  effect  could  be  virtually  eliminated  by 
applying  graphite  paint  around  the  edges  of  the  substrate. 


GRAPHITE 

PAINT 


Figure  4. 


APERTURE 


Experimental  Details  of 
Microstrip  Circuit  and 
Dig  Assembly 
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The  microstrip  configuration  described  presents 
considerable  advantages  over  the  open  resonator  which  was 
used  in  our  initial  experiments.  The  main  advantages  are 
outlined  belou: 

1.  For  similar  conditions  of  operation,  the 

signal  ue  obtained  is  approximately  lOdB’s 
higher  in  the  microstrip  case.  The  physical 
significance  of  this  fact  is  that:  though 


& 

t 


the  microstrip  resonator  has  a much  smaller 
Q (about  100  times  smaller),  the  microstrip 
resonator  is  very  much  better  at  confining 
tha  substrate  fields  to  a smaller  area  than 
the  open  resonator. 

2.  In  our  initial  experiments  with  the  open 

resonator,  drift  in  the  transmitter  frequency 
as  well  as  variations  in  the  cavity  resonator 
frequency  uas  the  major  source  of  signal 
instability.  This  effect  uas  largely  over- 
come by  the  incorporation  of  a frequency 
tracking  system  which  compensated  for 
fluctuations  in  either  the  cavity  resonant 
frequency  or  the  variations  in  the  transmitter 
output  frequency.  Owing  to  the  fact  that  the 
microstrip  resonator  is  an  inherently  low  Q 
device  as  compared  to  the  open  resonator,  this 
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pieces  a much  less  stringent  requirement  on 

the  tolerable  frequency  drift  due  to  the 

transmitter  or  the  cavity;  consequently 

the  frequency  tracking  system  was  not  required 

for  the  microstrip  case.  • 

3.  The  microstrip  resonator  is  not  only  substan- 

tially easier  to  fabricate  but  it  is  far 
simpler  to  alter  its  design  e.g.  the  coupling 
gaps  may  be  optimised  by  etching  off  the 
unwanted  copper.  In  addition,  the  micro- 
strip  version  is  considerably  lighter  and  more 
compact  than  the  open  resonator.  Microstrip 
technology  may  be  used  to  integrate  many  of 
the  components  used  in  our  system  onto  one 
single  substrate  (fig.  4.10):  this  is  not 

only  convenient  as  a space  saving  operation, 
but  has  the  added  advantage  that  the  resulting 
system  would  be  significantly  less  sensitive 

] 

to  both  thermal  and  microphonic  effects. 

4.  In  the  case  of  the  open  resonator,  air  coupling 
between  the  vibrating  object  and  the  lower 
resonator  plate  resulted  in  the  generation 

of  a slowly  varying  spurious  signal  as  tha 
object  was  translated  beneath  the  plav- 
The  effect  of  a fluctuating  spurious  signal 
is  to  reduce  the  overall  system  dynamic  range. 


Though  this  effect  could  be  eliminated  by 
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suitable  signal  processing  ue  would  loose 
some  of  our  main  signal.  This  problem  does 
not  occur  in  the  microstrip  case. 

5.  In  its  present  form,  the  microscope  is  limited 
to  the  inspection  of  flat  surfaces.  However, 
the  adoption  of  a microstrip  resonator  opens 
up  the  possibility  cf  inspecting  curved 
surfaces,  by  resorting  to  flexible  substrates. 

6.  One  significant  advantage  of  the  microstrip 
resonator  is  that  by  suitably  positioning  the 
aperture  in  the  ground  plane  ue  can  produce 
either  an  electric  or  a magnetic  dipole;  it 
would  be  very  much  harder  to  implement  this 
sort  of  change  in  the  open  resonator. 

4.2  Aperture  Design 

The  aperture  is  perhaps  the  single  most  important 
component  of  our  system,  Ue  have  already  demonstrated  that 

the  shape  and  size  of  the  aperture  determines  the  signal 
strength,  depth  of  field  and  the  resolution  ue  can  achieve. 

In  Section  2,3  ue  established  that  the  ability  of  the  micro- 
scope to  image  changes  in  the  properties  of  an  object  depended 
on  the  form  of  illumination  and  the  use  of  apertures  which 
were  small  as  compared  to  the  uavelength. 


43; 


The  significance  of  the  aperture  dimension  on  the 
resolution  has  already  been  established.  In  our  experiments 
we  have  used  circular  apertures  down  to  A/2GC  in  diameter 
and  have  demonstrated  a tuo  dimensional  resolution  capability 
of  better  than  A/375.  However,  if  ue  wished  to  improve  the 
resolution  by  a factor  of  say  10  by  decreasing  the  aperture 
diameter,  the  signal  would  be  reduced  by  a factor  of  30dB. 

A more  promising  approach  for  certain  NOT  applications,  such 
as  detecting  cracks,  which  usually  have  large  length  to  width 
aspect  ratios,  is  the  use  of  a rectangular  aperture.  To  a 
first  approximation  the  energy  stored  by  an  aperture  is 
proportional  to  its  area  (provided  the  significant  dimension 

is  very  much  less  than  the  wavelength).  If  we  keep  the 
areas  of  the  circular  and  the  rectangle  equal,  we  can  obtain 
roughly  the  same  signal.  By  appropriate  choice  of  the 
length  of  the  rectangle,  we  can  maintain  the  signal  level 
as  compared  to  the  circular  aperture  and  improve  the  resolu- 
tion by  decreasing  the  width.  In  addition,  the  time  :.'equired 
to  effect  a scan  would  be  reduced  by  the  factor  L/a,  where,  l, 
is  the  length  of  the  rectangle  and,  a,  is  the  aperture 
diameter.  The  penalty  for  improved  resolution  in  the  X 
direction  (Fig  4,11)  is  a corresponding  decrease  in  resolution 
in  the  Y direction.  Inspite  of  the  loss  of  resolution  in  one 
direction,  this  technique  seems  eminently  suitaole  for 
detecting  the  presence  of  cracks  in  metals. 


RE  SOLUTION(Y)  <C  L 
RESOLUTIONS)  «CW 


Figure  4.11 


Circular  and  Rectangular 
Aperture 
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5.  THE  EXPERIMENTAL  SYSTEM 

In  our  system  ue  are  using 

the  reflected  wave  from  the  object  and  consequently,  the 
detected  signal  consists  of  the  large  primary  reflection  from 
the  resonator  ground  plate  and  the  very  much  smaller  signal 
from  the  object.  By  only  detecting  the  signal  of  interest, 
ue  can  substantially  increase  the  input  pouer.  An  effective 
technique  to  discriminate  between  these  two  signals  is  to  make 
the  resonator  part  of  a bridge  circuit.  The  components  in 
the  other  arm  being  adjusted  so  as  to  suppress  the  large 
incident  reflection.  The  "contrast”  can  be  further  enhanced 
by  vibrating  the  object  at  an  acoustic  frequency  and 

subsequently  searching  for  this  component  in  the  receiver 
system. 


These  considerations  led  us  to  adopt  a suppressed- 
cai rier- single -sideband  arrangement;  the  main  elements  of 
the  system  which  are  shown  in  Fig.  5.1.  The  power  from  a 
10  GHz  impatt  oscillator  is  amplified  to  1 watt  by  a TWT 
amplifier  which  feeds  the  miciostrip  resonator;  the  atten- 
uator and  phase  shifter  in  the  bridge  circuit  are  adjusted 
to  suppress  the  carrier  so  as  not  to  overload  the  first  mixer 
The  single  sideband  modulator  (SSBM)  generates  a local 
oscillator  signal  which  is  locked  to  the  source,  and  which 
after  mixing  yields  an  i.f.  of  30  NH2.  This  is  then 
demodulated  using  a homodyne  technique  and  the  signal  at  Wffl 
is  finally  compared  in  a phase  sensitive  detectci  supplied 
with  a reference  from  the  acoustic  moaulator  source.  The 
output  of  the  phase  sensitive  detector  can  then  be  displayed. 


Figure  5.1  Experimental  Configuration  of  riicroacopa 


47. 


A theoretical  analysis  which  examined  the  noise  processes  associated  with 
transmitter  and  the  various  received  components  was  carried  out  in  order  to 
determine  the  major  sources  of  system  noise r 

The  nalysis  indicated  that  only  when  the  input  power  exceeded  100  watts 
would  one  be  dominated  by  a multiplicative  noise  process. — - 

Let  us  consider  the  "noise"  processes  associated 
with  the  transmitter.  We  have  established  that  the  effects 
of  transmitter  multiplicative  noise  may  be  neglected  in 
estimating  the  signal  to  noise  ratio*  We  have  verified  this 
experimentally  by  plotting  the  overall  system  signal  to  noise 
ratio  as  a function  of  input  power.  The  result  is  shown  in 
Fig.  5.4,  which  clearly  indicates  that  there  is  no  evidence  of 
transmitter  multiplicative  noise  for  powers  up  to  about  10  watts. 

We  have  found,  that  the  usable  dynamic  range  was 
substantially  lowered  on  account  of  slow  variations  of  the 
source  output  power  which  could  not  be  distinguished  from 
the  wanted  signal.  These  long  term  fluctuations  (occurring 
at  intervals  greater  than  1 second  in  our  case)  are  primarily 
due  to  thermal  effects  in  the  source  and  variations  in  the 
current  supplied  to  the  impatt  oscillator.  This  problem 
was  largely  overcome  by  the  inclusion  of  an  automatic  gain 
control  system  (AGC)  based  on  a PIN  modulator; 

Using  tnis  technique, 

we  have  reduced  fluctuations  in  the  output  power  from 
1:  5 x 10”  of  the  signal  to  about  1;  lo”5;  this  was  achieved 
for  variations  occurring  up  to  about  20  Hz.  Fig.  5.5  shows 
a trace  of  the  output  signal  from  the  system  with  and  without 
the  AGC  operating.  These  results  indicata  thot  even  with 
the  use  of  t..s  AGC  system,  we  are  limited  by  slow  variations 
in  the  transmitter  power.  The  actual  usable  dynamic  range 
is  therefore  only  50dB's  as  co,  -«d  to  the  short  te 
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Variations  of  the  object/aperture  spacing  while  translating  the  object 
results  in  the  generation  of  a spurious  signal  undistinguishable  from  the 
wanted  signal;  secondly,  any  deviation  from  parallelism  between  the  object 
surface  and  the  aperture  plane  results  in  a variation  of  the  output  signal 
which  limits  the  maximum  sensitivity  that  can  be  achieved. 

The  object  and  vibrator  were  scanned  below  the  .aperture  iut-A  X-Y  translation 
stage  (obtained  from  a travelling  microscope).  Even  slight  unevenness  in  the 

slide  mechanism  results  in  a spurious  signal;  for  example, 
a deviation  in  the  mean  separation  of  only  - 1 /J  would  cause 
a fluctuation  in  the  output  signal  of  about  - 0.1^.  It  is 
therefore  meaningful  to  distinguish  between  the  "static"  and 
the  "dynamic"  signal  to  noise  ratio.  The  static  measurement 
being  made  with  the  object  only  vibrating  under  the  aperture, 
and  the  dynamic  measurement  being  taken  while  the  object  is 
scanned  while  vibrating  beneath  the  aperture.  The  dynamic 
signal  to  noise  ratio  was  measured  for  various  positions  along 
the  translation  stage  while  the  specimen  was  being  scanned; 
the  final  usable  signal  to  noise  was  45  i 2dB. 

Let  us  consider  the  noise  added  by  the  receiver. 

It  was  found  that  the  first  mixer  was  the  major  source  of 
receiver  noise.  The  other  sources  of  noise  in  decreasing  * 
importance  were  the  second  mixer,  and  the  i.f,  amplifier. 

Fig,  5.6  shous  the  relative  noise  added  by  each  receiver 
component , 
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Figure  5,5 


Evaluating  Performance  of  AGC  - 
Signal  Output  uith  and  without 
AGC  in  Operation 
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Ue  have  nou  established  the  factors  which  are  limit- 
ing the  usable  dynamic  range  of  the  microscope.  These  are 
in  decreasing  importance:  variations  in  the  object/aperture 

separation,  fluctuations  in  the  source  power  (even  with  the 
AGC  in  operation),  receiver  additive  noise,  and  finally 
transmitter  multiplicative  noise. 

The  signal  to  noise  is  dependent  on  the  nature  of 
the  specimen,  the  aperture  dimension  and  the  final  system 
bandwidth.  The  resolution  and  detectability  is  also 
dependent  on  these  three  parameters.  The  overall  system 
performance  is  thus  necessarily  a compromise  between  the 
sensitivity,  resolution  and  the  time  required  to  effect  a 
scan.  The  following  example  illustrates  this  point. 

The  effects  of  fluctuations  in  the  transmitter  and 
variations  in  the  object/aperture  distance  can  be  reduced  by 
using  more  sophisticates  signal  processing  techniques  on  the 
records  obtained.  Consider  the  case  when  we  are  seeking  to 
detect  a specific  type  of  defect  such  as  a crack.  Ue  possess 
a priori  knowledge  that  the  cracks  are  small  and  hopefully 
rare.  This  information  can  be  used  selectively  to  filter 
out  very  high  fourier  components  that  are  physically  not 
possible,  as  well  as  filtering  out  the  lower  frequencies  which 
do  not  correspond  to  the  typical  crack  signature;  the  decon- 
volution procedure  described  in  Chapter  3 could  be  used  after 
this  preliminary  filtering  to  further  enhance  the  final  signal 
to  noise  ratio. 


6.  CHARACTERISING  THE  PERFORMANCE  OF  THE  MICROSCOPE 


In  this  chapter  the  results  obtained  from  a series  of  experiments  designed 
to  assess  the  performance  «f  the  microscope  and  correlate  these  results  with 
the  theory  developed  eai  jie"  are  presented. 

6.1.  Experimental  Operating  Conditions 

The  operating  conditions  under  which  the  measurements  were  made  are  as 
follows:  1 watt,  10  GHz  source  was  used  to  illuminate  the  iOO  u aperture  in  the 
microstrip  resonator-  the  samples  were  positioned  at  a mean  position 

of  about  50  y auay  from  the  microstrip  ground  plane,  and  were 

vibrated  (at  300  Hz)  at  an  amplitude  of  approximately  - 30  y . 

The  samples  were  translated  beneath  the  aperture  at  a constant 

speed  of  about  2.5  mm  per  minute.  In  all  the  experiments  the 

AGC  system  was  in  operation.  Tor  tests  on  metal  objects, 

the  aperture  uas  illuminated  so  as  to  produce  a magnetic  dipole. 

This  uas  because  the  signal  due  to  a magnetic  dipole  is  roughly 

6 d8's  higher  than  the  corresponding  signal  obtained  uith  an 

electric  dipole  for  metal  objects.  ___ _ . 

6.2  Signal  vs.  Distance  of  Plane 

Ue  have  developed  an  approximate  theory  of  the 
instrument  to  be  able  to  predict  hou  the  signal  varies  as  ue 
move  a conducting  plane  auay  from  the  aperture.  As  pointed 
out  earlier,  this  theory  becomes  progressively  more 

suspect  as  the  object  approaches  uithin  one  aperture  diameter. 
Nevertheless,  ue  find  some  encouraging  agreement  betueen  the 
theory  and  experiment  as  shoun  in  Fig.  6,1.  The  best  fit 
here  is  obtained  using  tuo  images  in  the  range  0.25  to  1.0 
aperture  diameters  auay  from  the  diaphragm.  It  should  be 
noted  that  the  theory  uas  derived  for  a square  aperture  uhereas 
in  practice  ue  are  using  a circular  one;  and  that  the 
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theoretical  and  experimental  results  were  matched  at  a 
distance  of  one  aperture. 

6.3  Signal  vs.  Aperture  Diameter 

In  Chapter  3 ue  established  that  the  signal  uas  pro- 
portional to  the  cube  of  the  aperture  diameter.  Experiments 
were  performed  to  test  the  validity  of  this  theory.  Apertures 
of  350,  300,  250  and  200  microns  uere  fabricated  in  the  ground 
plane  of  the  microstrip  resonator  and  the  signal  obtained  as 
a metallic  object  uas  vibrated  at  a constant  amplitude  under 
each  aperture  uas  recorded. 

The  test  object  uas  vibrated  at  a fixed  amplitude 
of  about  - 30  ^j,  uith  a mean  separation  of  approximately  50  y 
auay  from  the  aperture  plane.  In  order  to  compare  the  signal 
obtained  from  each  aperture,  ue  must  normalise  these  results. 


The 

normalised  signal  to  noise  ratio  for  these  four  apertures 
is  plotted  in  Fig,  6.2.  The  theoretical  and  practical 
results  agree  uithin  the  limitations  of  the  experiments,  and 
the  calculated  correcting  factors. 
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6.4  Besolution  Tests 

The  attainable  resolution  is  limited  by  the  aperture 
dimension,  the  dynamic  range  of  the  detection  system,  and  the 
properties  of  the  sample  ue  are  investigating.  Ue  have  carried 
out  tests  on  metallic  gratings  consisting  of  a 0.2  y aluminium 
film  on  glass  substrates  to  assess  the  resolutions  ue  can 
achieve  for  such  an  "ideal"  object.  Fig.  6.3(a)  shows  a 
scan  across  a grating  of  line-width  125  y using  an  aperture 
diameter  of  300  jj.  Fig.  6.3(b)  shows  a similar  scan  obtained 
with  an  aperture  diameter  of  200  y , for  a grating  having  a 
line-width  of  only  80  jj;  this  result  demonstrates  a resolution 
capability  of  better  than  ^/375  with  the  present  instrument. 
Experiments  on  polished  glass  substrates  are  admittedly  at  least 
one  step  away  from  realistic  n.d.t.  applications;  houever, 
there  are  some  problems  which  do  involve  the  detection  of 
microcracks  on  smooth  metallic  surfaces;  and  the  results 
obtained  provide  a useful  guide  to  assess  system  performance. 

6.5  Results  on  51ots  in  Metal  Objects 

Ue  have  developed  an  approximate  theory  to  be  able 
to  predict  how  the  signal  varies  as  a function  of  the  depth 
and  width  of  a slot  located  symmetrically  below  a square 
aperture.  Quantitive  experiments  were  performed  using  a 
circular  aperture  to  test  the  validity  of  this  theory. 

Variable  depth  spark  machined  slots  in  brass  of  the 
same  width  were  investigated.  The  results  are  shown  in 
Fig.  6.4(a),  where  we  have  plotted  the  maximum  signal  from 
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tuo  slots  of  different  cross-sections  against  the  depth  of  the 
slot,  together  with  the  theoretically  predicted  results. 

The  theory  and  experiment  suggest  that  ue  are  not  able  to 
provide  information  of  the  depth  of  the  slot  if  the  UID1H/ 
DEPTH<£  1.  However,  for  certain  n.d.t.  applications  it  is 
only  necessary  to  indicate  the  presence  of  the  defect  and 
for  these  cases  the  instrument  could  be  used  to  advantage, 

Experiments  were  also  carried  out  on  slots  of 
variable  width  and  constant  depth.  The  results  are  shown 
in  Fig.  6.4(b)  and  are  compared  with  the  curve  computed  from 
the  theory  developed  m Section  3.4.  While  the  theory 
predicts  that  the  signal  should  be  proportional  to  the  square 
of  the  width,  the  experimental  results  indicate  that  the 
dependence  on  width  is  in  fact  more  nearly  linear.  The  gap 
between  the  theory  and  experiment  has  been  previously  explained; 
any  improvement  in  the  theory  would  entail  performing  a detailed 
quasi-static  field  calculation. 


Figurs  6.4(a)  Signal  Strength  vs.  Depth  of  Slot  for 
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7.  APPLICATION  TO  NON-DESTRUCTIVE  TESTING 


In  this  chapter  ue  assess  the  potential  application 
of  our  microscope  to  the  non-destructive  evaluation  of  surface 
topography  and  surface  properties  of  materials. 

Microuave  methods  are  naturally  the  most  direct  for 
testing  materials  such  as  dielectric^  ^ ^ and  ferrite 

materials^ used  in  microuave  components.  Their  use  is 
also  increasing  in  other  n.d.t.  problems  such  as  the  evaluation 
of  bonded  films^  of  surface  texture^*  \ and  detection  of 
fatigue  cracks  in  metals. 


Most  of  these  methods  rely  on  the  mode  conversion  or  scattering  due  to 
the  defect  on  the  illuminated  surface.  Such  methods  can  be  made  remarkably 
sensitive  in  detecting  single  defects.  One  could  also  imagine  that  they  might 
come  into  use  for  the  statistical  evaluation  of  surface  finish.  However, 
compared  to  other  non-destructive  testing  techniques  such  as  X-rays,  acoustics 
and  luminescent  paints,  they  suffer  from  the  poor  ability  to  locate  the  precise 
position  of  one  of  a series  of  defects.  This  lack  of  locateability  and 


resolution  derives  from  the  wavelength  of  the  microwaves. 
There  is  thus  a potential  need  for  a technique,  which  gives 
resolution  capability  which  is  very  much  greater  than  with 


a 
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other  methods  at  the  same  frequency.  In  order  to  ascertain 
the  potential  areas  in  which  our  microscope  would  be  particularly 
advantageous,  we  have  carried  out  tests  on  a variety  of  samples 
which  may  be  classified  into  five  basic  categories: 

1*  The  detection  of  micro-cracks  in  solid 
metal  objects, 

2.  The  detection  of  hairline  breaks  in  thin 
metallic  films  on  dielectric  substrates. 

3.  Recording  changes  in  the  relative  per- 
meability of  ferrite  samples. 

4.  Imaging  small  variations  in  the  electric 
permittivity  of  dielectric  objects. 

5.  Characterising  the  properties  of  thin 
films. 

In  this  chapter,  ue  shall  present  the  results  obtained  on 
each  of  the  five  classes  of  defects  outlined  above  and  discuss 
the  possible  applications  of  the  microscope  to  realistic  n.d.t. 
situations. 


I 
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Detection  of  Micro-Cracks  in  Metals 


Initial  experiments  were  directed  towards  examining 
the  potential  use  of  the  instrument  in  the  evaluation  of 
surface  topography.  Fig.  7.1  shows  a typical  result  of  a 
defect  in  a sample  of  steel,  where  a scan  of  the  flaw  is 
compared  with  a cross-section  of  the  object  as  viewed  through 
a metallurgical  microscope.  It  is  interesting  to  see  that 
some  of  the  irregularities  in  the  observed  record  do  in  fact 
correspond  to  the  actual  shape  of  the  object;  though  naturally 
extended  and  distorted  due  to  the  fact  that  the  record  is 
essentially  a convolution  of  the  defect  with  the  aperture 
and  electronic  system  functions.  This  implies  that,  provided 
the  defects  ue  wish  to  image  are  larger  than  the  resolution 
limit  of  the  instrument,  the  microscope  could  be  used  in  the 
statistical  examination  of  surface  topography  if  suitable 
calibration  techniques  were  used  to  relate  the  signal  to 
actual  height  variations. 


The  experiments  on  fabricated  slots  in  metals  presented 
in  Section  b.5  suggested  that  the  instrument  was  capable  of 
detecting  defects  which  were  significantly  smaller  than  the 
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aperture  dimension.  Ue  have  carried  out  experiments  to 
enable  us  to  estimate  the  minimum  detectable  defect.  Pig. 

7.2  shows  a trace  from  a 2 y wide,  3 mm  deep  steel  slot, 
together  with  a scanning  electron  micrograph  of  the  slot; 
which  was  fabricated  by  joining  two  highly  polished  steel 
pieces.  In  this  case,  the  sample  was  translated  at  a 
constant  speed  of  0.1  mm/10  mt. 

Despite  the  precautions  taken,  we  could  not  be  sure 
that  ue  might  not  be  detecting  some  characteristic  which  would 
not  be  present  in  a true  crack.  For  this  reason,  samples 
uere  obtained  by  the  use  of  a fatigue  machine  on  steel  speci- 
mens. The  smallest  defect  ue  have  been  able  to  detect  was  a 
fatigue  crack  believed  to  be  less  than  1 y wide  in  a tool 
steel  sample.  The  result  is  shown  in  Fig.  7.3,  which  also 
shows  a scanning  electron  micrograph  of  the  crack.  Compared 
to  other  microwave  crack  detectors^  ^ >^,13  these  results 
are  more  than  two  orders  of  magnitude  better  in  terms  of  both 
sensitivity  of  detection  as  well  as  in  resolution  capability. 
It  should  be  pointed  out  that  though  ue  can  detect  cracks 
which  are  significantly  smaller  than  the  aperture  diameter, 
ue  would  not  be  able  to  resolve  two  1 y cracks  separated  by  a 
distance  less  than  the  current  resolution  limit,  80  y in  our 
case  (using  a 200  y aperture). 

In  some  n.d.t.  applications,  the  metal  surface  to  be 
inspected  is  coated  with  a protective  layer,  which  is  often 
essentially  a dielectric  coating  e.g.  a randome.  It  was 
therefore  of  interest  to  see  whether  the  microscope  could 
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detect  a crack  beneath  such  a protective  layer.  Fig.  7.4 
shows  that  though  the  crack  produced  signal  is  diminished  due 
to  the  50  jj  dielectric  coating,  the  flaw  is  still  readily 
detected.  This  result  uas  obtained  with  an  aperture  of 
diameter,  700  jj. 

7.2  Detection  of  Breaks  in  Thin  P.etallic  Films 

Ue  have  demonstrated  the  ability  of  the  microscope 
to  resolve  thin  metallic  gratings  having  line  widths  less  than 
80  u,  .A  closely  related  situation,  and  one  uhich  comes  very 
much  closer  to  a class  of  material  testing  problems,  is  to 
locate  a single  hairline  crack  in  an  otherwise  continuous 
metal  film.  Fig.  7.5  shows  a result  for  a deliberate  micro- 
crack  of  uidth  2 jj,  uhich,  it  is  seen,  is  still  within  the 
detection  limit  of  the  instrument.  Similar  results  wire 
obtained  for  cracks  in  thin  metallic  films  on  silicon  sub- 
strates. 


In  order  to  be  useful  as  a n.d.t.  instrument,  the 
microscope  should  be  capable  of  discriminating  between  two 
surface  defects:  first,  a break  in  a thin  metallic  film  on 

a dielectric  substrate  and  secondly,  a surface  scratch  of 
similar  dimensions.  Fig.  7.6  shows  that  though  the  scratch 
(or  surface  depression)  is  not  detected  by  our  microscope, 
the  bleak  is  clearly  resolved.  In  this  experiment  we  chose 
the  aperture  fields  to  correspond  to  a purely  magnetic  dipole 
for  which  the  near-field  stored  electric  energy  is  very  small, 
so  that  the  si'vai  in  the  case  of  the  break  in  the  metallic 
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film  arises  predominantly  from  gradients  in  the  distribution. 

In  the  case  of  the  surface  scratch,  currents  are  set  up  in  the 
depression  and  consequently,  the  resulting  field  is  not  altered 
significantly  to  effect  a detectable  change  in  the  signal. 

The  ability  to  detect  hairline  cracks  in  thin  metallic  layers 
is  already  on  the  edge  of  direct  usefulness  in  applications 
such  as  thin  film  circuitry, 

Ue  have  performed  quantitive  experiments  to  be  able 

to  determine  hou  the  signal  varies  as  a function  of  the  uidth 
« 

of  the  slit  in  a thin  metallic  film  on  a dielectric  substrate. 
Variable  uidth  slits  were  fabricated  in  a 0.2  y aluminium 
film  on  a glass  substrate.  Fig.  7.7  shows  that  the  signal 
is  very  nearly  a linear  function  of  the  slit  uidth.  Ue  have 

& 
t 

made  no  attempt  to  predict  this  relationship  rigorously  as 
this  would  require  a very  complex  field  calculation. 


Nore  theoretical  work  is  requireo  to  predict  the  smallest 
break  in  a thin  metallic  film  on  a dielectric  substrate. 


Variations  in  Dielectric  Con stant  and 


magnetic  Permeability 

Besides  its  use  in  high  resolution  inspection  of  surface  topography, 
the  microscope  is  capable  of  providing  information  about  the 
electric  and  magnetic  properties  of  the  test  specimen.  it 
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uas  shoun  that  if  an  electric  dipole  uas  used,  the  microscope 
would  be  very  much  more  sensitive  to  changes  in  the  relative 
dielectric  constant  than  to  changes  in  the  magnetic  perme- 
ability of  the  object.  Similarly,  for  a magnetic  dipole, 
the  instrument  would  be  more  sensitive  to  variations  in  the 
permeability  of  the  sample.  This  provides  the  basis  of  a 
novel  technique  to  selectively  image  variations  in  either 
the  dielectric  constant  or  the  magnetic  permeability  of  an 
object  with  a very  high  resolution.  In  order  to  corroborate 
this  theory,  ue  have  carried  out  a series  of  experiments  on 
dielectrics,  ferrites  and  finally  cast  iron  alloys. 

The  microscope  uas  used  to  record  spatial  variations 
of  the  relative  dielectric  constant  uith  an  electric  dipole. 
Fig.  7.9  shows  a scan  from  the  edge  of  a polystyrene  sample 
to  a perspex  sample,  where  it  can  be  seen  that  a variation  of 
only  2.5$  in  relative  dielectric  constant  uas  easily  detected. 
The  figure  suggests  that  with  the  present  instrument, 
variations  of  only  0.3$  in  dielectric  constant  could  be 
resolved.  This  makes  the  microscope  potentially  useful  in 
the  high  resolution  testing  of  dielectric  substrates  such  as 
microstrip  substrates. 


It  is  important  to  appreciate  that  since  the  depth 
of  field  is  itself  very  much  smaller  than  the  aperture  diameter, 
that  we  can  "see"  surface  changes  which  uould  be  extraordinar- 
ily difficult  to  detect  in  any  other  way  at  the  same  frequency 
(simply  because  for  handlable  specimens,  the  "filling  factor" 


i 
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would  be  very  poor).  This  opens  up  quite  a new  form  of 
microwave  measurements  - to  reveal  the  behaviour  of  dielectric 
and  ferrite  thin  films,  Ue  have  carried  out  tests  on  thin 
dielectric  films  to  verify  this  important  property.  M45 
Dielectric  films  serve  many  essential  fi  notions  in  modern 
solid  state  device  technology 

it  is  important  to  measure  their  quality  and  to  determine 
the  presence  and  nature  of  any  imperfections  that  may  give 
rise  to  device  failure  so  that  appropriate  corrections  can 
be  made  in  thar  processing. 

(13  ) 

Recently,  Giallorenzi  et  cl'  ' and,  Pitt  end 

( 14  ) 

Planku'  , have  demonstrated  the  possibility  of  using  ion- 
exchange  in  glasses  to  fabricate  low  loss  optical  waveguides. 
The  depth  of  the  guide  and  the  diffused  ion  concentration 
profile  are  important  parameters  in  determining  the  perform- 
ance of  the  guide;  however,  they  are  difficult  measurements 
to  make  in  practice  accurately.  Ue  have  carried  out 
preliminary  tests  on  a sample  of  soda-glass  substrate  with 
an  ion-exchanged  (silver  in  our  case)  film  grating  which  was 
only  about  5 p deep.  The  results  are  shown  in  fig.  7.10; 
where  it  is  seen  that  the  thin  dielectric  film  is  readily 
detected.  The  change  in  dielectric  constant  of  the  film  to 
the  substrate  was  estimated  to  be  of  the  order  of  10$. 

Recently,  the  suitability  of  elect rostrictive 
materials,  such  as  potassium  tantalate  niobate  (KTM)  have 


( 15  ) 

been  investigated  ' as  substrate  materials  for  the  genera- 
tion of  surface  acoustic  waves.  An  important  parameter  in 
determining  the  efficiency  of  generation  is  the  dependence 
of  dielectric  permittivity  on  the  applied  bias  voltage.  Ue 
are  not  aware  o;  any  published  information  on  the  dielectric 
behaviour  of  KTN  as  a function  of  electric  fxeld;  it  was 
therefore  of  interest  to  see  whether  the  microscope  could 
record  changes  in  the  dielectric  constant  on  application  of 
a bias  field.  Using  the  configuration  shown  in  fig.  7.11(a) 
ue  recorded  the  signal  as  a function  of  the  distance  from  the 
two  electrodes.  The  asymmetric  behaviour  of  the  signal  with 
position  on  the  top  surface  (Fig.  7,ll(a)  ) indicates  that 
dielectric  permittivity  is  dependent  on  the  polarization  of 
the  electric  field  with  respect  to  the  crystal  orientation, 
fig.  7.11(b)  shows  how  the  dielectric  constant  varies  with 
applied  bias  when  the  aperture  was  positioned  midway  between 
the  two  electrodes. 

These  results  are  potentially  useful  in  assessing  the  perform 
ancc  of  the  material  for  possible  applications  in  SAU  genera- 
tion, Tne  results  must  be  interpreted  carefully  if 

they  are  to  provide  useful  data:  as  for  large 

perturbations  in  the  dielectric  constant,  the  signal  is  no 
longer  a linear  function  of  the  dielectric  permittivity. 

The  possibility  of  seeing  magnetic  contrast  was 
also  explorer).  By  using  an  aperture  Distribution  correspond 
ing  to  a magnetic  dipole,  we  have  shown  the  "visibility"  of 
a ferrite  sample,  fig,  7.12  shout  a scan  across  a metal/ 
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alumina/spinel  ferrite  substrate.  The  dielectric/ferrite 
boundary  is  clearly  resolved.  Having  demonstrated  the 
ability  of  "seeing"  magnetic  contrast,  it  was  of  interest 
to  assess  the  sensitivity  of  the  instrument  to  changes  in 
the  magnetic  permeability  of  a ferrite  sample.  Preliminary 
experiments  uere  performed  on  a lcw-f requency  tuning  coil 
ferrite  to  measure  the  change  in  magnetic  permeability  with 
applied  magnetic  field.  The  ferrite  uas  poled  in  one  direction 
by  applying  a large  d.c.  field  and  the  change  in  permeability 
uith  an  increasing  field  in  the  opposite  sense  uas  recorded. 

The  result  is  shown  in  fig.  7.13(a);  on  integration  of  this 
curve,  we  obtain  the  B vs.  H for  the  ferrite  at  microwave 
frequencies.  Though  we  would  sxpect  some  variation  of  the 
microwave  permeability  (|jf)  as  a function  of  the  magnetic 
field,  it  is  not  clear  what  the  effective  ^rvs.  H curve  should 
look  like.  However,  the  real  value  of  the  experiment  is  in 
that,  it  clearly  demonstrates  that  by  using  a magnetic  dipole 
aperture,  the  microscope  can  image  small  variations  in  the 
magnetic  permeability.  Ue  believe  that  the  microscope  has 
direct  application  in  imaging  small  variations  in  magnetic 
thin  film  materials  such  as  used  in  bubble  memories^  ’ ' 

to  assess  for  example  the  material  uniformity  which  is  an 
important  parameter  in  determining  the  performance  of  the 
ferrite. 

One  further  application  of  the  microscope  may  be 
in  the  field  of  metallurgy.  Conventional  optical  metailo- 
graphic  techniques  involve  three  basic  steps:  first,  prepara- 

tion of  a section  surface;  secondly,  the  development  of  a 
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structure  usually  by  a chemical  etching  process;  and  finally, 
the  actual  observation  and  recording  of  the  structure*  This 
is  both  destructive  and  time  consuming.  Provided  the  resolu- 
tion of  our  microscope  is  adequate  for  resolving  the  details 
ue  wish  to  investigate,  it  may  be  used  to  advantage  to 
characterise  certain  alloys.  Ue  have  investigated  three 

commercially  used  cast  iron  alloys.  These  alloys  uere  of 

( T8) 

interest  as  they  contained  secondary  phase  particles'  , 
the  electric  and  magnetic  properties  of  which  varied  over  the 
specimen.  Fig.  14(a)  and  (b)  show  the  results  of  scans  across 
a sample  of  hyper-eutectic  white  cast  iron  and  a specimen  of 
phosphoric  white  cast  iron  respectively,  together  with  a 
scanning  electron  micrograph  of  the  area  under  inspection. 

In  Fig.  14(a),  the  most  prominent  dark  cementite  needles 
(Fe^C)  are  resolved  as  is  the  small  nest  cf  flake  graphite. 

In  Fig.  14(b),  the  large  30  u surface  crack  is  easily  detected 
and  recorded  as  a dip  in  the  signal,  whereas  the  two  dark 
pearlite  (ferrite  and  cementite)  regions  show  up  as  an 
increase  in  the  relative  signal. 

Ue  have  demonstrated  that  besides  its  possible  use 
in  the  non-destructive  testing  of  metal  surfaces  to  detect 
sub-micron  cracks,  the  microscope  has  potential  application 
in  providing  information  about  the  surface  properties  of  both 
ferrite  and  dielectric  thin  films  which  would  be  very  difficult 
to  detect  in  any  other  way.  However,  before  the  instrument 
can  be  used  in  practical  n.d.t.  si;uations,  more  work  in  cal- 
ibrating the  sensitivity  of  the  microscope  to  specific 
perturbations  would  be  necessary. 
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B.  SURFACE  DEFECT  DETECTION  USING  ONE-DIMENSIONAL  ACOUSTIC  HOLOGRAPHY 
1.  INTRODUCTION 

The  basis  of  the  method  is  to  illuminate  the  surface  to  be  tested  with  acoustic 
surface  waves,  launched  at  one  end  of  the  sample  and  then  to  record  the  amplitude 
and  phase  distribution  along  a scanning  line  at  the  other  end  of  the  sample. 

This  record  can  be  regarded  as  a hologram  which  can  subsequently  be  reconstructed 
so  as  to  reveal  features  at  any  point  on  the  surface.  At  the  time  of  the 
original  proposal  which  led  to  the  current  work,  we  had  succeeded  in  making  such 
a hologram  with  a phase  sensitive  laser  probe,  and  reconstructing  the  record  using 
conventional  laser-optic  techniques.  We  had  already  recognized  though  that  there 
would  be  very  substantial  advantages  in  using  a computer  reconstruction  of  the  data. 
Since  the  record  is  one-dimensional  the  reconstruction  computation  is  not  arduous  - it 
takes  only  a few  seconds  on  an  IBM  360  65.  We  have  also  succeeded  in  effecting  the 
reconstruction  on  quite  a small  laboratory  computer. 

The  work  which  has  been  accomplished  on  this  project  has  to  a large  extent 
been  published  - at  least  in  preliminary  form.  Since  the  published  statement  is 
suitably  concise,  we  are  reprinting  the  three  relevant  papers  in  the  following. 

The  first  paper  describes  the  work  immediately  preceding  the  award  of  the  contract 
and  the  initial  work  on  computer  reconstruction  carried  out  during  the  first 
quarter  of  the  work.  The  most  significant  results  relevant  to  n.d.t.  appear 
in  figure  7,  which  shows  the  successful  computer  reconstruction  of  a number  of  small 
gold  dot  "defects",  and  the  first  successful  determination  of  a slowness  character- 
istic in  figure  11.  Ir  addition  it  gives  a concise  description  of  the  laser  probe 
system,  and  also  of  the  effects  of  noise  on  the  reconstruction  of  images. 

Subsequently  a considerable  amount  of  further  effort  was  devoted  to  improving 
the  laser  probe  both  with  regard  to  its  electrical  as  its  mechanical  performance. 

This  led  eventually  to  the  much  superior  results  on  slowness  determination,  described 
in  the  second  publication,  section  3. 

In  section  4 we  present  a letter  publication  which  reports  a major  improvement 
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in  the  holographic  defect  detection  technique.  It  rests  on  the  linearity  of  the 
process  which  permits  one  to  subtract  the  record  of  a known  good  sample  from  that 
obtained  with  a sample  which  may  have  a defect.  It  led  to  the  identification  of 
a gold  dot  “defect"  only  30  pin  diameter  and  0.4  p thick  a defect  which  was  totally 
buried  in  the  “noise"  without  the  use  of  this  subtraction  technique.  We  believe 
that  this  technique  will  prove  of  major  importance  in  the  development  of  these 
methods. 

In  section  5 we  briefly  summarize  the  results  obtained,  and  attempt  an 
assessment  of  their  potential  importance  to  n.d.t.  applications. 
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Optical  laser  probes  have  for  long  been  used  for  measurements  on  acoustic  surface  wave  devices. 
At  frequencies  below  a few  hundred  megahertz  one  can  use  a probing  system  in  which  the  focused 
spot  is  small  as  compared  with  the  acoustic  wavelength.  This  enables  one  to  derive  fine-grained  ampli- 
tude and  phase  information  on  the  surface  wave  field.  Some  examples  of  the  diagnostic  use  of  a par- 
ticular instrument  of  this  design  will  be  presented. 

Since  the  probe  gives  complete  amplitude  and  phase  information,  a scan  transverse  to  an  acoustic 
surface  wave  beam  contains  in  principle  a complete  diffraction  limited  record  of  the  two-dimensional 
field,  right  up  to  the  transmitting  transducer.  We  can  regard  the  line  scan  as  a one-dimensional 
hologram  of  the  two-dimensional  object  space.  The  hologram  can  be  reconstructed  either  optically 
or  on  a computer.  Experiments  relating  to  both  methods  will  be  presented.  The  technique  is  seen 
as  having  potential  for  the  nondestructive  testing  of  surfaces.  It  is  not  necessarily  confined  to  plane 
surfaces. 

Effective  reconstruction  on  the  computer  implies  a detailed  knowledge  of  the  slowness  surface 
of  the  substrate.  Conversely  one  ca”  derive  the  slowness  surface  from  a comparison  of  the  complex 
distribution  along  two  line  scans.  The  technique  therefore  represents  a method  for  measuring 
the  slowness  surface,  and  is  capable  of  attaining  very  high  accuracies.  Preliminary  results  of  some 
measurements  will  be  presented. 


I.  iNTRODUCTION 

Acoustic  surface  waves  are  readily  launched  and  received;  it  is  less  easy  to  observe 
their  propagation  from  transmitter  to  receiver.  The  displacement  amplitudes  involved 
are  small -typically  several  orders  of  magnitude  less  than  one  angstrom.  Nevertheless, 
a number  of  successful  probe  detection  systems  have  been  devised.  In  one  class  of 
device,  one  uses  a stylus  directly  attached  to  a piezoelectric  transducer  [ 1 ] ; when  probing 
nonmetallized  substrates  which  aie  themselves  piezoelectric,  a simple  metal  stylus  picks 
up  a voltage  which  can  be  directly  amplified  [2] . Such  probes  have  the  inherent  advan- 
tage of  simplicity. 

A second,  class  is  based  on  the  use  of  lasers,  in  which  the  phase  modulation  oc- 
casioned by  the  surface  displacement  is  in  some  way  detected.  Laser  probes  can  them- 
selves be  divided  into  a number  of  categones  depending  on  whether  the  interaction  can 
be  regarded  as  a diffraction  process,  extended  over  many  acoustic  wavelengths  f3] , or  as 
a phase  modulation,  confined  to  at  most  one  half  acoustic  wavelength.  The  present  work 
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is  based  on  the  use  of  the  latter  type,  in  a form  first  described  by  Whitman  and  Korpel 
[4]-  Laser  probes  have  the  immediate  advantage  of  causing  a negligible  perturbation  of 
the  acoustic  surface  wave,  and  of  avoiding  any  possible  scratching  of  the  surface.  More- 
over they  work  as  well  on  insulating  as  on  conducting  surfaces.  The  price  is  a greater 
degree  of  complexity. 

The  laser  probe  as  originally  described  was  designed  for  amplitude  measurements 
only  [5] . For  many  of  our  purposes  we  were  particularly  interested  in  the  phase  distri- 
bution. For  example,  in  the  far  field  of  a transducer  the  relevant  information  is  almost 
entirely  phase  information.  We  therefore  decided  to  make  our  probe  capable  of  mea- 
suring phase  as  well  as  amplitude.  The  basic  design  adopted  has  been  previously  de- 
scribed [6] , and  analyzed.  We  will  in  Section  II  give  a very  brief  description  of  the 
apparatus  in  its  present  improved  form. 

In  this  paper  we  want  to  address  ourselves  to  the  use  of  laser  probes  for  two  different 
classes  of  application.  The  first  is  concerned  very  directly  with  the  development  of 
acoustic  surface  wave  devices,  and  measurements  to  diagnose  sources  of  malfunction.  In 
Section  III  we  will  give  one  example  of  the  use  of  the  laser  probe  for  such  purposes.  The 
probe  is  seen  to  be  of  particular  importance  in  connection  with  devices  which  are  truly 
two  dimensional  in  character  (i.e.,  where  the  propagation  effects  in  the  transverse  direc- 
tion are  of  essential  significance  to  the  basic  working  of  the  device). 

The  laser  probe  can,  however,  also  be  used  in  a somewhat  different  manner  for  the 
characterization  of  a complete  two-dimensional  field  in  a single  scan.  Just  as  a two- 
dimensional  holographic  record  can  store  information  about  a three-dimensional  object 
space,  a one-dimensional  record-a  single  laser  scan -can  contain  complete  information 
(within  the  usual  resolution  limits)  on  a two-dimensional  field.  The  process  by  which 
one  can  reconstruct  the  scene  from  the  scan  information  has  been  discussed  by  a number 
of  workers  in  the  context  of  ordinary  (i.e.,  three-dimensional)  optics,  and  acoustics. 

In  the  case  of  surface  waves,  it  leads  to  the  reconstruction  of  the  two-dimensional  fields, 
and  in  particular  to  the  detection  of  scattering  sources  within  that  field  by  observation  at 
a distance.  The  concept  is  therefore  one  which  lends  itself,  in  principle,  to  the  non- 
destructive testing  of  components.  Some  preliminary  experiments  to  this  end  are 
described  in  Section  IV. 

In  the  case  of  anisotropic  materials  the  reconstruction  of  the  field  from  a scan  re- 
quires an  accurate  knowledge  of  the  slowness  curve  of  the  substrate.  Where  this  is  not 
known  with  sufficient  accuracy,  the  method  suggests  a new  measurement  technique:  it 
involves  the  recording  of  two  scans  separated  by  a distance.  The  relationship  between 
the  spectra  of  the  two  scans  then  leads  directly  to  the  delineation  of  the  slowness  curve. 
A discussion  of  the  method,  and  some  preliminary  results  are  also  presented  in  Section 
IV. 


II.  LASER  PROBE 

The  basic  design  of  the  probe  [6]  is  shown  in  Figure  1 . It  can  be  regarded  as  a modi- 
fied Michelson  interferometer  in  which  the  beam  splitter  is  replaced  by  a Bragg  cell 
thereby  providing  a carrier  frequency  for  the  signal.  The  introduction  of  such  a carrier 
frequency  brings  some  advantages  in  the  subsequent  signal  processing;  above  all  it  ensures 
that  spurious  optical  path  length  variations  arc  not  detiiment.il  provided  that  their  spec- 
trum does  not  fall  into  the  signal  processing  bandwidth,  centered  on  tens  of  megahertz. 
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Fig.  1.  Laser  heterodyne  probe. 


The  system  is  therefore  remarkably  free  from  disturbances  which  would  otherwise  be 
occasioned  by  microphonics,  small  temperature  gradients,  etc. 

In  our  first  experiments  to  obtain  phase  discrimination  [7] , we  vibrated  the  reference 
. mirror  at  the  acoustic  signal  frequency  fs,  while  retaining  simple  square  law  detection  of 
the  signal.  This  technique  allows  one  to  use  a bridge  method,  where  the  phase  of  the 
measured  signal  is  derived  from  the  phase  change  setting  in  the  reference  arm  to  obtain 
a null  in  the  detection  system  [8] . It  is  a method  which  works  effectively  but  which  is 
difficult  to  extend  so  as  to  give  a continuous  phase  record  as  the  sample  is  scanned. 

The  natural  alternative  is  to  use  a fixed  reference  mirror,  but  with  a phase-coherent 
cetection  system.  The  frequency  component  or  interest  after  photodetection  is 
(2/b  - /,),  where  fD  is  tee  Bragg  cell  frequency.  This  suggests  that  one  could  derive 
a local  oscillator  signal  by  doubling  the  Bragg  cell  drive  frequency.  However,  it  turns 
out  that  by  so  doing,  one  looses  the  insensitivity  to  path  length  fluctuations  at  low  fre- 
quencies which  is  such  a vital  factor  in  the  realization  of  a stable  instrument.  Fortunately 
there  is  another  source  of  2 fy  in  one  component  of  the  photodetected  laser  current. 

This  component  contains  the  same  optical  phase  information  as  the  ( lfB  - fs ) com- 
ponent. After  mixing  one  has  therefore  once  again  restored  the  immunity  to  optical 
phase  jitter.  .These  considerations  then  led  to  the  detection  scheme  shown  in 
Figure  2. 

We  have  found  that  even  in  this  situation  one  can  experience  signal  fluctuations 
arising  from  hydrodynamic  streaming  instabilities  in  the  Bragg  cell.  With  careful  design, 
and  control  of  temperature  of  operation  this  effect  can  be  largely  avoided. 

The  present  instrument  is  designed  to  work  at  a series  of  bands  up  to  100  MHz.  No 
major  obstacles  to  its  extension  to  frequencies  two  or  three  times  greater  are  envisaged. 

In  operation  the  sample  is  mechanically  scanned  under  the  laser  beam;  the  vector 
voltmeter  output  can  be  displayed  directly  on  an  XY  recorder.  Alternatively,  and  this 
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Fig.  2.  Coherent  laser  probe  detection  system. 


the  signal  can  be  fed  directly  to  a minicomputer  which  digitizes  and  in  some  cases  pro- 
cesses the  signal. 

It  is  of  particular  relevance  briefly  to  examine  the  phase  resolution  of  the  instrument. 
At  a frequency  of  60  MHz,  the  acoustic  wavelength  in  quartz  and  in  LiNb03  is  of  the 
order  of  50  jx.  The  laser  probe  can  be  focused  to  less  than  5 n-  The  phase  detection  sys- 
tem has  a resolution  limit  which  is  better  than  10°.  Since  one  can  expect  to  obtain  a 
resolution  which  is  less  than  the  phase  width  of  the  spot,  it  is  this  latter  figure  which  is 
probably  the  limiting  factor.  We  can  therefore  achieve  a resolution  at  this  frequency  of 
around  A/20.  With  a sample  length  of  5 cm,  we  have  a total  path  length  of  103  wave- 
lengths. The  instrument  should  therefore  be  capable  of  measuring  velocity  with  an 
accuracy  of  around  5 X 1CT5 . The  performance  could  be  improved  if  longer  samples  are 
available.  In  practice,  the  attainment  of  such  accuracies  also  implies'considerable  effort 
on  removing  or  discriminating  against  bulk  wave  and  other  spurious  signals. 


\ 

j 
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III.  DEVICE  DIAGNOSTICS 

Measurements  of  amplitude  distribution  along  acoustic  surface  wave  devices  have 
been  published  by,  many  authors.  In  some  cases  these  have  served  to  reveal  aspects  of 
performance  which  it  would  have  been  hard  to  explore  in  any  other  way  [9] . The  basis 
and  many  examples  of  such  measurements  are  presented  in  a recent  review  article  [3] . 

In  this  section  we  will  give  a single  recent  [ lOj  example  of  the  use  of  the  laser  probe  In 
a situation  where  both  the  amplitude  and  phase  distribution  was  of  importance. 

The  problem  is  concerned  with  transduction  onto  a topographic  waveguide.  The 
difficulty  of  direct  transduction  [11,  12]  stems  from  the  fact  that  the  “aperture”  of  the 
guide  is  less  than  one  wavelength,  so  that  the  transducer  has  to  be  fabricated  at  the  very 
edge  of  the  guide,  a process  which  it  is  difficult  to  accomplish  using  normal  photolitho- 
graphic techniques.  We  have  sought  to  effect  this  coupling  using  a conventional  inter- 
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Fig.  3.  (a)  Thin  film  to  wedge  guide  directional  coupler.  The  amplitude  distribution,  recorded  at 
60  MHz  shows  the  progressive  transfer  of  energy  to  the  wedge  guided  wave,  (b)  Phase  plot  along  the 
edge  of  the  wedge  guide.  There  axe  39  wavelengths  in  the  2 mm  scan  corresponding  to  a velocity 
of  3088  m/sec. 


OPTICAL  AND  ACOUSTICAL  MICRO  ELECTRONICS 


87. 


to  the  topographic  waveguide.  Figure  3 shows  the  system  used  in  a particular  experiment 
to  couple  60  MHz  waves  onto  a 64'  angle  wedge  waveguide.  The  gold  waveguide  has  a 
thickness  so  that  its  velocity  is  depressed  10%  below  the  Rayleigh  velocity  for  the  Ycut 
‘ sample  (in  the  * propagation  direction).  Computations  subsequently  showed  [13J  that 
the  wedge  wave  for  this  angle  is  only  2%  below  the  Rayleigh  velocity ; however  since  in 
j this  case  the  thin  Film  guide  actually  touches  the  edge  of  the  wedge  guide,  the  coupling 
is  very  strong;  a substantial  transfer  of  energy  tiierefore  takes  place,  in  spite  of  the  lack 
of  synchronism.  The  probe  amplitude  distributions.  Fig.  3(a),  clearly  show  the  manner 
in  which  the  guide  wave  is  launched,  and  the  decrease  in  the  amplitude  of  the  thin 
film  guide  mode. 

Since  the  synchronous  condition  is  central  to  the  operation  of  this  device,  it  is  clearly 
important  to  have  an  accurate  knowledge  of  the  wedge  guide  velocity.  This  is  obtained 
using  the  laser  probe,  simply  by  recording  the  phase  as  the  probe  is  scanned  along  the 
edge  of  the  wedge  guide.  The  phase  record  is  shown  in  Figure  3(b).  Although  in  this 
case  only  a relatively  short  s an  was  used,  the  measured  velocity  is  determined  within  an 
accuracy  of  a few  parts  in  103 , in  good  agreement  with  the  value  computed  by 
LaGasse  f 13] . 


IV.  ONE-DIMENSIONAL  HOLOGRAPHY 

We  will  now  turn  to  the  use  of  the  laser  probe  for  the  recording  of  one-dimensional 
holograms.  It  is  possible  to  record  such  a hologram  as  a one-dimensional  record  on  a 
photographic  plate,  and  then  subsequently  to  effect  the  reconstruction  optically.  An 
example  of  this  technique  [14]  is  shown  in  Figure  4(a).  The  “object”  consisting  of  two 
gold  dots  on  a quartz  substrate  was  illuminated  with  surface  waves  at  60  MHz,  and  laser 
probed  at  a distance  of  400  wavelengths.  The  scan  was  recorded  in  the  form  of  a binary 
phase  hologram  which  is  shown  in  Fig.  4(b);  the  optical  reconstruction  which  clearly 
resolves  the  two  dots  is  shown  in  Figure  4(c). 

The  optical  reconstruction  however  represents  a troublesome  additional  step.  More- 
over it  leads  to  the  distortion  arising  from  the  difference  of  lateral  and  longitudinal 
magnifications  in  the  ratio  of  the  acoustic  to  the  optical  wavelength  [15].  It  is  therefore 
natural  to  seek  to  effect  the  reconstruction  numerically.  The  main  difficulty  encountered 
in  doing  this  in  the  case  of  two-dimensional  holograms  - the  sheer  magnitude  of  the 
numerical  task-is  not  a significant  factor  in  the  case  of  a one-dimensional  record.  More- 
over one  can  use  a reconstruction  wavelength  identical  to  that  used  in  the  recording,  so 
that  object  distortion  is  totally  avoided.  There  is  in  addition  the  opportunity  to  use  a 
number  of  techniques  which  can  be  adopted  to  make  the  fullest  possible  use  of  any  a 
prjori  knowledge  of  the  “object”  which  is  sought.  Finally  it  is  possible  to  use  the  tech- 
nique for  anisotropic  media,  even  in  the  case  where  the  slowness  curve  cannot  be  ade- 
quately represented  by  a parabola. 

The  process  with  which  we  are  concerned  is  equivalent  to  “inverse  diffraction” 

| [ 1 6-20] , and  there  are  a number  of  publications  which  have  specifically  applied  the 

technique  to  the  reconstruction  of  normal  (three-dimensional)  acoustic  holograms 
[21-24] . Our  case  is  somewhat  simpler  in  that  we  are  concerned  with  only  two  dimen- 
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Fig.  4.  (a)  Sample  used  in  two-dimensional  holography,  (b)  Binary  phase  hologram  recorded  at 
60  MHz,  400  wavelengths  from  object,  (c)  Optical  reconstruction  of  object. 
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j A.  Reconstruction  Theory 

i 

j Figure  5 shows  an  “object”  distribution  U0  (x)  at  2 = 0 and  an  “image”  distribu- 

| tion  (//  (x)  at  z = L.  A plane  surface  wave,  propagating  in  the  0 direction  will  have  a 
j spatial  dependence  of  exp  (-  j (ax  + 0z)),  where 

a = k (0)sin  (0) , 

P = k (0)  cos  (0) , (1) 

and  k (0)  is  the  wavenumber.  We  can  obtain  an  expansion  of  the  object  distribution  in 
terms  of  the  plane  wave  spectrum, 

C am 

(/»M=  /o  (a)  e~,ax  da . (2) 

The  waves  corresponding  to  values  of  a greater  than  the  maximum  value  attained  by 
k (0)  sin  (0)  arc  evanescent,  and  can  certainly  be  excluded  from  our  calculations,  since 
we  will  be  concerned  with  object-image  distances  of  many  wavelengths.  We  have  there- 
fore used  the  limits  a;„  and  a/„  corresponding  to  the  maximum  values  of  a for  0 positive 
and  negative,  respectively.  Provided  that  the  magnitude  of/0  (a)  drops  to  low  values 
by  the  time  wc  approach  these  limits,  we  may  regard  Kq.  (2)  as  a simple  Fourier  trans- 
form, with  the  inverse 

fo(a)'-  f U0(x)cinxdx.  (3) 

J-00 

Since /o  (a)  can  be  interpreted  as  the  amplitude  of  the  plane  wave  corresponding  to  0,  we 
can  [25]  immediately  write  down  the  corresponding  spectral  density  at  the  image  plane 

fi  (a)  - Jo  (a)  e . (4) 

This  is  a standard  procedure  in  isotropic  diffraction  theory.  It  applies  equally  in  the 
anisotropic  case,  provided  that  0 (a)  is  a single  valued  function  of  a.  Wc  will  make 
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this  assumption  for  thcmoment,  but  revert  to  the  problem  raised  by  double  valued 
functions.  We  can  then  at  once  take  the  transform  of  //  (a)  and  hence  obtain  the  image 
distribution 


l//(x)  = f Jo  (ot)e  ML+ax)  tja . 


(5) 


The  procedure  can  be  reversed  so  as  to  yield  the  “object”  distribution,  from  a known 
“image”  distribution.  Using  precisely  the  same  method  one  obtains 


i 

c ain 

1 

Uo  (x)  = 1 or Jf  1 

UI(x)e>axdx\da., 

(6) 

I 

J-oo  j 

In  the  case  of  some  materials,  the  function  0 (a)  is  not  single  valued  for  the  whole 
range  of  1 0 1 < nil.  An  example  is  presented  in  Fig.  6 which  shows  a (0)  as  a function 
of  0 for  Z-cut  rutile,  0 = 0 corresponding  to  the  Y axis.  We  see  that  for  0 greater  than 
43°,  there  arc  two  and  for  some  ranges  of  a , three  possible  values  of  0,  and  hence  through 
Eq.  (1),  two  possible  values  for  (3 ; there  is  no  longer  a unique  result  for  the  integral  of 
Equation  (6).  The  physical  significance  of  this  conclusion  is  that  the  distribution  U0  (x) 
is  simply  not  uniquely  defined  by  Ui  (x). 

It  is  interesting  to  relate  this  conclusion  to  a simple  case,  where  the  object  distribu- 
tion is  due  to  two  plane  waves,  having  the  same  value  of  a,  but  two  different  values  of 
0,0!,  and  02.  lf/li  and  A2  are  the  corresponding  complex  amplitudes 


and 


U0(x)  = (. A , +A2)e’l**t 


Uj(x)  = (/l,  +A2  L) e~'ax 


(7) 


(8) 


Clearly  the  relationship  between  A{  and  A2  cannot  be  derived  from  a measurement 
of  only  U0  or  only  Uf.  However  if  we  measure  both  U0  and  UI  we  have  two  complex 
equations  for  the  two  complex  amplitudes,  and  we  can  proceed  to  solve  for  both  A x and 
A2.  In  fact  we  see  that  if  we  measure  any  two  distributions  displaced  by  a known  length 
in  the  z direction,  we  can  obtain  the  complete  plane  wave  spectrum,  and  hence  derive 
the  distribution  at  any  other  plane;  We  see  therefore  that  in  order  to  reconstruct  an 
object  in  this  situation,  we  need  to  measure  two  image  distributions  displaced  by  an 
arbitrary  distance  in  the  z direction.  Finally  it  is  interesting  to  note  that  the  spectral 
intensity  of  the  two  distributions  of  Eqs.  (7)  and  (8),  are  no  longer  identical,  as  they  in- 
evitably are  when  0 (a)  is  single  valued. 

The  experimental  work  so  far  completed  relates  to  the  single  valued  case,  to  which 
we  will  now  return.  The  formulation  of  Eq.  (f>)  is  still  a step  or  two  removed  from  the 
reality  of  an  experiment,  where  Uj  (x)  is  obtained  as  a sampled  rather  than  as  a con- 
tinuous function.  The  function  is  sampled  at  N equally  spaced  points,  displaced  by 
Ax.  It  is  convenient,  particularly  when  using  Fast  Fourier  Transform  techniques,  to 
compute  the  spectrum  also  at  N points,  displaced  by  Au.  Finally  the  object  distribution 
is  again  reconstructed  in  terms  of  N samples  with  the  same  spatial  increment  Ax.  The 
discrete  Fourier  transform  version  of  Eq.  ((>)  then  reads 

</v/2)-t  -to  , * . <n/2> 

U0(q  Ax)  = £ eHPpr.-pqAa*x>  y U{  (^'Ajc)  clpq  aaa x , (9) 

d=-(NI2)  q 


£ Uj  (r/’Ax)  elpq  AaAx  , 
(A >12) 
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where 

Pp  = {k  (O)2  - (pAa)2  )1{2  . 

If  we  use  the  Fast  Fourier  Transform,  we  must  restrict  N to  a oinary  number.  Also  we 
must  adopt  a specific  relationship  between  Aa  and  Ax 

2ir 

' AaAx  = — • (10) 

n 

i 

With  Eq.  (10),  Eq.  (9)  can  be  written  in  the  final  form 

U0(qAx)  = (yV£  * ^ ‘ U,  {q  Ax)e*2vp<*'IN)  . (11) 

P-~(N  ID  q'=-(H/2) 

| This  is  the  form  in  which  the  computations  were  carried  out. 

In  the  case  of  such  a one-dimensional  record  the  use  of  Fast  Fourier  Transform 
methods  is  not  mandatory,  and  indeed  it  would  be  entirely  possible  to  use  diffraction 
formulations  other  than  plane  wave  expansions.  It  is  however,  for  our  purpose,  the  par- 
ticular merit  of  the  plane  wave  expansion  that  the  anisotropy  of  the  medium  is  so  readily 
included,  at  least  as  long  as  0 (a)  remains  single  valued.  One  can  see  at  once  that  in  ordei 
to  achieve  an  accurate  reconstruction  of  the  object  distribution  the  form  of  the  aniso- 
tropy, k (0)  must  be  known  with  considerable  precision,  as  the  transform  contains  a 
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phase  factor  of  the  order  of  A'  (0)  L , and  L may  in  practice  be  a large  number  of 
wavelengths. 

B.  Reconstruction  of  Experimental  Data 

Preliminary  experiments  have  been  carried  out  using  a quart/  substrate  and  using  a 
distribution  of  gold  dots  to  represent  the  “object."  Figu?„  7(a)  shows  the  arrangement 
and  the  position  of  the  various  scan  records.  The  gold  dots  had  diameters  of  ^ mm,  and 
a thickness  of  8()00  A.  Hie  “object"  distribution  St , taken  immediately  adjacent  to  the 
line  of  dots  is  shown  in  Fig.  7(b),  and  the  “image”  distribution  .V2  in  Figure  7(c).  The 
computer  reconstruction  in  accordance  with  Hq.  ( 1 1 ) is  shown  first  for  a line  S3  displaced 
by  5 mm  from  the  object  (Fig.  7(d))  and  finally  focussed  onto  the  object  line  itself 
Figure  7(c).  It  is  seen  that  the  three  dots  arc  well  resolved. 

In  nondestructive  testing  applications  (Section  IV.D)  it  may  prove  valuable  to  be 
able  to  reconstruct  unknown  defects,  in  the  presence  of  known  “defects”  (e.g.,  the  edge 
of  the  sample).  Since  all  the  manipulations  involved  in  the  reconstruction  process  are 
linear,  it  is  clearly  possible  to  add  and  subtract  records.  Experiments  to  demonstrate 
this  possibility  are  in  progress. 

The  experiments  arc  in  an  early  phase,  and  the  limits  of  the  detection  sensitivity  of 
the  laser  probe  has  yet  to  be  explored.  However  in  the  next  section  we  will  discuss  some 
preliminary  considerations  which  bear  on  this  issue. 


C.  Detection  and  Resolution  Sensitivity 

The  resolution  and  detection  sensitivity  arc  limited  both  by  the  normal  limits  of 
resolution  of  any  optical  instrument,  as  by  the  finite  accuracy  and  signal-to-noise  ratio 
achieved  in  the  laser  probe.  However  it  is  of  interest  to  note  that  in  principle  the  detec- 
tion limit  may  be  far  below  the  Rayleigh  criterion,  provided  that  one  ha:  some  a priori 
knowledge  of  the  nature  of  the  object.  It  is  a common  place  that  one  can  see  stars  which 
one  cannot  resolve.  Very  simple  considerations  indicate  that  this  situation  is  one  which 
one  might  well  encounter  in  practice.  Suppose  that  the  sample  is  illuminated  by  a total 
power  P,  applied  over  a width  IV.  The  power  density  is  then  P/W,  and,  if  we  are  probing 
in  the  near  field,  will  be  of  the  same  order  in  the  image  scan  position.  Suppose  that  we 
have  a totally  absorbing  object,  width  6 IV,  spaced  D from  the  image  line.  This  object 
can  then  be  represented  as  giving  rise  to  a cylindrical  wave  having  a (negative)  power  of 
Ph  1V/1V.  The  amplitude  ratio  between  this  scattered  wave  and  the  largely  unperturbed 
illumination  is  then  R,  where 

R2  = (5  1V/2t r/7) . 


The  scattered  wave  will  at  different  points  have  a phase  of  ± rr/2  relative  to  the  main 
wave,  leading  to  total  phase  variations  of  the  resultant  50,  where 

50  - (25IV/7rD),/2  . (12) 


Suppose  that  we  are  able  to  measure  0 with  an  accuracy  of  1°,  and  that  D = 100X.  We 
then  find  that  5 IV  ~ A/20.  The  resolution  must  of  course  always  be  of  the  order  of  a 
wavelength.  However,  if  one  is  looking  for  a sparse  set  of  defects,  or  even  for  the  possible 
presence  of  a single  one,  we  see  that  it  may  prove  possible  to  detect  very  small 
disturbances. 

In  order  to  obtain  a more  detailed  picture  of  the  probe  performance  required  to 
achieve  a particular  resolution,  we  have  simulated  the  random  errors  in  a theoretical 
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reconstruction  of  three  point  sources,  spaced  by  four  wavelengths.  We  have  computed 
(for  an  isotropic  substrate)  the  image  distribution  at  a plane  100  wavelengths  from  the 
object.  Using  Eq.  (1 1)  with  N = 512,  we  have  then  reconstructed  the  distribution  at 
various  positions  back  to  the  original  object  plane  (Fig.  8(a),  with  the  actual  distribution 
at  the  object  plane  shown  separately  in  Figure  8(b)).  The  point  source  distribution  has 
been  broadened,  a i expected  from  the  finite  spatial  bandwidth  of  the  assumed  system. 

We  can  now  simulate  tin  inadequacies  of  the  laser  probe  by  adding  random  signals 
to  the  theoretical  sampled  probe  output.  In  conducting  such  a simulation  we  must 
choo$c  the  standard  deviation  for  amplitude  and  phase,  which  we  wish  to  impose.  At 
ti  e same  time,  we  must  simulate  the  approximate  spatial  frequency  of  the  assumed 
noise.  Tire  highest  frequency  noise  is  obtained  if  we  add  random  uncorrelated  signals 
to  each  adjacent  cell.  We  can  simulate  lower  spatial  frequency  disturbances  by  choosing 
a random  signal  for  every  nth  sampling  interval,  and  making  a linear  interpolation  for 
tlie  intermediate  sampling  points. 

Figures  9(a)  and  (b)  show  the  effect  of  adding  pure  amplitude  noise,  with  no  corre- 
lation between  adjacent  sampling  points  (i.e.,  with  n = 1 ).  It  is  seen  that  even  for 
AA/A,  the  standard  deviation  referred  to  the  peak  value,  as  large  as  1.2  the  three  sources 
are  still  just  resolved.  Figures  9(c)  and  (d)  siiow  similar  results,  again  for  //  = 1 , for  pure 
phase  noise.  Here  we  find  that  the  reconstruction  still  allows  for  identification  of  the 
sources  for  A 0 as  large  as  80°. 

The  addition  of  uncorrelatcd  noise  at  adjacent  sampling  steps  provides  a good  model 
of  the  deleterious  effects  arising  from  errors  in  the  analog-to-digilal  converter,  and  also 
for  the  effects  of  relatively  high  frequency  noise  in  the  detection  circuitry.  However 
just  because  there  is  no  correlation  between  adjacent  steps,  which,  in  this  case  arc 
separated  by  about  A/5,  the  noise  spectrum  lies  mainly  in  the  evanescent  part  of  the 
spectrum,  and  in  the  propagating  spectrum  corresponding  to  waves  which  travel  at  large 
angics  to  the  normal  so  that  their  effect  on  the  reconstruction  is  minimized.  Noise 
with  the  same  standard  deviation,  but  correlated  over  several  sample  points  (i.e.,  for 
n > I),  has  a larger  deleterious  effect  on  the  quality  of  reconstruction.  This  is  apparent 
from  Fig.  9(c),  where  we  have  used  the  ratio  of  the  peak  in  the  reconstruction  of  the 
point  sources  to  the  amplitude  of  the  largest  competing  spurious  maximum  as  a some- 
what arbitrary  measure  of  the  damage  done  by  the  presence  of  the  phase  noise.  The 
figure  shows  how  this  ratio  varies  with  the  value  of//.  The  dashed  line  indicates  the 
ratio  for  the  reconstruction  in  the  absence  of  noise. 


D.  Nondestructive  Testing 

Experiments  on  a polished  quartz  surface  arc  admittedly  at  least  one  step  away  from 
a realistic  n.d.t.  situation.  However  there  are  at  least  some  problems  which  do  involve 


the  detection  of  microcracks  in  smooth  surfaces,  though  more  generally  the  techniques 
with  which  we  are  concerned  are  usable  on  components  which  are  not  optically  fiat; 
we  have  cairied  out  a wide  range  of  meusuiements  on  metal  samples  using  an  earlier, 
version  of  the  piesent  laser  probe  [7], 

The  reconstruction  of  the  gold  dots  (Fig.  7(e))  is  a primitive  example  of  the  “defect" 
detection  possibilities.  It  is  clear  howevei  tha'  one  could  improve  the  detection  sensi- 
tivity by  resorting  to  more  complex  detection  ol  small  delects.  In  many  cases  theielore 


the  task  will  be  concerned  with  the  identification  ol  point  sources.  This  a priori  knowl- 
edge of  the  natuie  of  the  “object”  should  greatly  facilitate  in  both  its  detection  and 
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location.  At  present  we  arc  digitizing  the  data  on  a minicomputer,  but  carrying  out  the 
reconstruction  on  a larger  machine.  In  fact  it  is  entirely  feasible  to  carry  out  the  whole 
process,  in  (almost)  real  time.  Now  the  assumed  point  defect  will  give  rise  to  a phase  dis- 
tribution having  the  character  of  a Fresnel  /.one  plate.  If  then  we  correlate  the  scan  data 
with  such  a phase  distribution,  we  will  for  each  choice  of  range  obtain  a correlation  peak, 
the  largest  occurring  for  the  range  corresponding  to  the  actual  object  location.  The  pro- 
cessing gain  associated  with  this  technique  can  be  quite  large,  and  could  lead  to  a 
detection  sensitivity  which  is  substantially  better  than  that  indicated  by  the  simple  cal- 
culation leading  to  Equation  (12). 

Having  suspected  the  presence  of  a defect  at  a particular  location,  one  could  then 
seek  confirmation  by  the  use  of  a second  scan,  somewha'  neater  to  the  defect.  Indeed 
it  might  prove  convenient  always  to  work  with  more  than  one  scan,  and  directly  to 
process  the  results  obtained  from  each. 

In  practice  few  of  the  objects  that  one  would  wish  to  test  would  be  planar.  Provided 
that  the  shape  of  the  object  is  such  that  it  can  be  illuminated  by  surface  waves,  and  pro- 
vided that  at  one  end  there  is  a suitable  location  for  a scanning  line,  the  laser  scan  tech- 
nique is  in  principle  applicable.  Ia  such  a case  one  could  not  usually  be  able  to  effect 
direct  computer  reconstructions.  However,  it  mtght  be  possible  to  map  the  surface  by 
measuring  the  “image”  scan  distribution  when  particular  points  on  the  sample  surface 
arc  excited.  In  effect  one  would  be  obtaining  an  experimental  Green’s  function  for  the 
particular  sample.  This  could  then  be  compared  with  potentially  defective  samples,  to 
identify  and  locate  faults. 

E.  Slowness  Curve  Determination 

The  reconstructions  for  the  quartz  samples  discussed  in  Section  IV.  B,  involved  an 
accurate  knowledge  of  the  slowness.  In  fact  earlier  attempts  where  the  slowness  data 
was  inaccurate  to  the  extent  of  a few  parts  in  l O'1  led  to  reconstructions  marred  by  vio- 
lent ripples.  This  suggests  that  the  technique  offers  a means  for  the  accurate  determina- 
tion of  slowness  curves.  In  principle,  the  procedure  is  quite  simple.  One  measures 
Uq  (x)  and  £//  (x)  and  hence  derives  /0  (a)  and  J)  (a).  From  Eq.  (4)  we  then  have 


j3  /.  = cos 


[.//(a)  , //(ft) 
\ 2 L /o(«)  /,)  (a) 


+ 2nN 


- F(a)  + 2t jN  , 


where  F( a)  is  defined  by  Eq.  ( 1 3),  and  N is  an  integer.  Now 

a = k (0)  sin  (0) , 


P = --  ( F(a ) + 2nN)  = k (0)  cos  (0) , 

so  that 


[A'(0)J2  = cv’  t-  'j  {F\ cvH  2nNf  , 


/'(«)  + 2ni\ 
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Fig.  10.  Amplitude  scans  for  slowness  determination,  (a)  Immediately  beyond  random  array  of  dots. 

(b)  100  \ beyond  random  array. 


Equations  (14)  and  (15)  arc  applicable  to  the  case  when  0 (a)  is  single  valued.  In  order 
to  obtain  usable  results  it  is  naturally  essential  that  the  object  spcct/um  contains  com- 
ponents over  the  whole  range  of  0 for  which  the  slowness  surface  is  required.  A simple 
solution  to  this  requirement  is  to  use  a very  narrow  aperture  transducer;  however  this  is 
at  the  expense  of  the  total  available  signal.  One  can  also  usi  a large  suitably  curved  trans- 
ducer. An  alternative  which  we  have  adopted  in  our  prebn, inary  experiments  is  to  use 
a random  array  of  small  gold  viots  illuminated  by  a normal,  wide  interdigital  transducer. 

Figure  10(a)  shows  the  amplitude  scan  as  measured  just  beyond  the  second  (of  two) 
rows  of  dots,  and  Fig.  10(b),  the  amplitude  at  the  image  plane.  We  then  computed  k ( 0 ) 
as  a function  of  0 from  Equations  ( 14)  and  (15).  The  results  arc  shown  in  Fig.  1 1 , and 
arc  compared  with  the  curves  computed  from  elastic  constant  data  by  Siobodnik  [26] . 
The  form  of  the  slowness  characteristic  is  reproduced,  though  the  accuracy  is  modest. 

It  is  believed  that  the  gap  between  the  achieved  accuracy  and  that  which  we  believe  to 
be  possible  with  the  apparatus  is  primarily  due  to  the  fact  that  we  had  taken  only  minimal 
precautions  against  bulk  and  other  spurious  waves.  Futute  experiments  are  being 
designed  with  this  problem  very  much  in  mind. 

In  principle  the  accuracy  which  one  could  achieve  ought  to  approach  the  basic 
j^putitant^ftreeisioti  discussed  in  Section  II.  One  can  also  hope  to  improve  the  basic 
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Fig.  11.  Velocity  results  for  T-cut  quartz  at  0°  to  the  x direction  compared  with  Slobodnik  et  al  [26] . 


accuracy  by  making  post-recording  corrections.  These  can  be  based  on  one  known  rela- 
tionship between  the  object  and  image  tecords  the  fact  that  for  single  valued  0 (a),  the 
amplitude  of  the  spectral  density  functions, /(a)  must  be  identical.  As  a simple  example, 
consider  the  case  of  a scan  recorded  subject  to  a slow  phase  drift  in  the  recording  system. 
This  means  that  instead  of  recording  the  true  distribution  U (x)  one  will  be  recording 
U'  (x),  where 


(/'(*)=  U(x)CikrX  , 
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ft- 
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It 
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where  kr  indicates  the  magnitude  of  the  phase  drift.  This  will  lead  to  a transform /'  (a), 
related  to  the  true  transform  /(a)  by 


/'  (a)  = /(<*  + kr) 


06) 


Tlius  we  see  that  uniform  drift  leads  to  a simple  shift  in  the  spectrum.  We  can  readily 
shift  the  two  spectra, /0  (a)  and  //  (a),  to  bring  them  into  optimum  coincidence,  before 
computing /(a)  through  Eq.  (4)  and  then  k ( 0 ) through  Equation  (14).  Removal  of  a 
linear  drift  does  not  exhaust  the  possibilities  of  correction  for  more  complex  distortions; 
one  can  anticipate  that  further  improvements  in  the  recorded  data  can  be  obtained  in 
this  way. 

The  question  m mains  as  to  how  one  can  extend  this  technique  for  the  case  of  an 
anisotropic  substrate,  for  which  0 (a)  is  no  longer  single  valued.  The  procedure  is  in  this 
case  somewhat  more  involved  but  not  in  any  basic  way  more  complicated.  Essentially 
one  needs  to  rccoid  additional  scans,  the  number  depending  on  the  multiplicity  of  the 
0 (a)  function.  For  a double  valued  function  we  record  one  additional  scan  (i.e. , three 
in  all)  and  additional  scans  for  each  unit  increase  in  the  multiplicity.  This  then  provides 
enough  information  to  allow  a unique  determination  of  k {0). 
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3.  SURFACE  ACOUSTIC  WAVE  SLOWNESS  SURFACE  MEASUREMENT 
H.K.  Wickramasinghe  and  E.A.  Ash 

ABSTRACT.  Surface  wave  velocity  curves  can  be  measured  with  high  precision  using  a phase  sensitive 
laser  probe  The  method  does  not  imply  any  assumptions  about  the  wave  field  distribution  used.  Accuracies 
of  a few  parts  in  10*  have  been  attained,  and  suggest  that  the  technique  may  also  prove  valuable  in  the 
determination  of  fundamental  elastic  constants. 


Introduction 

Critical  design  problems  involving  surface  acoustic 
waves  require  an  accurate  knowledge  of  the  slowness 
curve  of  the  substrate.  For  some  commonly  used  materials 
such  as  Lithium  Niobate  and  Bizmuth  Germanium  Oxide, 
these  surfaces  are  not  known  with  great  accuracy.  In 
this  paper,  we. report  a technique  whereby  the  slowness 
surface  can  be  measured  with  high  accuracy  over  a wide 
angular  range.. 

The  technique  involves  the  measurement  of  two  trans- 
verse complex  field  distributions  spaced  by  a known 
distance  using  a phase  sensitive  laser  probe('). 

Relating  the  spectra  of  these  two  scans  leads  directly 
to  the  slowness  curve  of  the  substrate.  The  evaluation 
makes  no  assumptions  about  the  distribution  of  the 
incident  surface  acoustic  wave  field. 

A detailed  analysis  of  the  theory  involved  was 
presented  in  a previous  publication!2).  Following  a 
brief  outline  of  this  theory,  we  will  present  experi- 
mental results  demonstrating  accuracies  of  a few  parts 
in  1C)4  for  some  commonly  used  materials.  Finally  we 
will  discuss  the  capabilities  of  the  technique  in  the 
light  of  our  results. 


l Igure  1 shows  as  "object"  distribution  U0(x)  at 
Z = 0 and  an ‘image"  distribution  tii(x)  at  Z=L.  A plane 
surface  wave  propagating  in  the  e direction  will  have  a 
spatial  dependance  of  exp-j(ax+Bz)  where 

a = K(e)  sine 

e = K(e)  cose  (l) 

and  K(0)  is  the  wave  number.  We  can  expand  U„(x)  in 
terms  of  the  plane  wave  spectrum 


UQ(x)  = / - VQ>e'jaX  do 
_am 


The  waves  corresponding  to  values  of  a greater  than  the 
maximum  value  attained  by  K(o)  sine  are  evanescent, 
and  can  certainly  be  excluded  from  our  calculations 
since  we  are  concerned  with  propagation  distances  of 
many  wavelengths.  We  have  therefore  used  the  limits 

a+  and  a*  corresponding  to  the  maximum  values  of  a for 
in  m 

6 positive  and  negative  respectively.  Provided  that 
the  magnitude  of  f0{a)  drops  to  low  values  by  the  time 
we  approach  these  limits,  we  may  regard  equation  (2) 
as  a simple  Fourier  transform  with  the  inverse 


f0M  ■ f y*)ejax  d* 


since  f0(a)  can  be  interpreted  as  the  amplitude  of  the 
plane  wave  corresponding  to  e,  we  can ( 3 ) immediately 
write  down  the  corresponding  spectral  density  at  the 
"image"  plane, 

f|W  = f0(«)e‘jBL  (45 


This  is  a standard  procedure  in  isotropic  diffr- 
action theory.  It  applies  equally  to  the  anisotropic 
case  provided  e(a)  is  a single  valued  function  of  a(h 
We  can  then  at  once  take  the  transform  of  fi(a)  and 
hence  obtain  the  image  distribution, 

Uj(x)  » /“m  f0(o)e*i(BL+ax}da  (5) 

‘"i 

To  calculate  an  accurate  image  distribution  U.(x) 
from  a measured  object  distribution  U0(x),  the  form  of 
the  slowness  curve  K(o)  must  be  known  with  considerabl. 
precision  since  the  integrand  in  equation  (5)  contains 
a phase  factor  BL  and  L may  in  practice  be  several 
hundred  wavelengths.  Clearly,  one  can  reverse  this 
process  to  obtain  an  accurate  representation  of  the 
slowness  curve  by  measuring  Uj { x) ; the  accuracy 
achieved  being  directly  proportional  to  L. 

To  do  this,  one  met  ures  U0(x)  and  U.(x)  and  henc- 
derives  f fa)  and  f,(o).  From  equation  (4)  we  then 


BL  = cos’1  J 


>;(«)  t 

_V8>  +fM 


= F(q)  + 2nN 

where  F(a)  is  defined  by  equation  (6)  and  N is  an 
integer.  One  then  finds  that  « 


K(e)2  « a2  + £[F(a)  + 2nNj  2 

9 „ tan'1  f gL  : -1 
[F(a)+2»NJ 


The  theory  as  it  stands  does  not  enable  us  to  max 
an  absolute  measurement  of  the  slowness  surface,  since 
the  value  of  N is  not  known.  The  value  of  N can  be 
determined  however  if  we  know  K(o)  for  any  arbitrary 
value  of  o.  Tms  information  could  be  obtained  experi 
mentally  by  making  a separate  measurement  in  th» 
particular  direction  of  interest.  For  the  sake  of 
convenience  we  measure  K(0).  The  value  of  N follows 
immediately  from  equation  (7). 

2nN  = K(0)  L - F(0)  (9) 

The  velocity  determination  effected  by  a measure- 
-ment  of  tne  number  of  wavelengths  in  a known  distance 
inplies  far  field  conditions.  Using  the  wider  trans- 
ducers this  condition  is  not  met,  and  can  lead  to  pnas 
errors  which  are  significant  in  the  context  of  the  hi. 
accuracies  attainable  by  the  method.  One  can  avoid  *• 
need  for  such,  a separate  measurement  by  resorting  to 
a third  scan  Uj(x)  at  Z = L ‘ . From  equation  (51  we 
have 
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Uj(x) 


e’joX  da 


were  chosen  to  be  the  same  as  those  for  the  previous 
experiment.  The  experiment  enabled  us  to  obtain  the 
(10)  slowness  curve  over  a wider  (3-?°)  angular  range. 

Discussion  and  Conclusion 


substituting  f„r  $ froiu  (a)  gives 
+ 

Uj(x)  = /“m  f#(.)e-^F<“)^  e'jaX  da  (11) 


Ir.  order  to  determine  the  integer,  N,  we  then  simply 
vary  N until  we  obtain  the  best  correlation  betweer  the 
measured  distribution  Uj(x)  and  that  calculated  from 
equation  (11). 

Experimental  Results 

It  is  essential  that  the  object  spectrum  contains 
components  over  the  whole  range  of  : for  which  the 
slowness  surface  is  required.  A simple  solution  to 
this  requirement  i*.  to  use  a very  narrow  aperture 
transducer.  However,  this  is  at  the  expense  of  th-* 
total  available  signal.  An  alternative  is  ..o  use  a 
random  array  of  small  gold  dots  illuminated  by  a wid. 
interdigital  transducer.  The  latter  alternative  has  an 
advantage  in  that  one  can  subtract  the  effect  of  bulk 
and  other  spurious  waves  by  recording  Up( x ) and  Uj ( x) 
with  and  without  the  random  gold  dots(  3). 

Figure  2 shows  the  measured  slowness  curve  for  YX 
Quartz  with  a 600  A overlay  of  chrome,  compared  with 
the  theoretical  curve  for  non-metall ized  Quartz  com- 
puted by  Slobodnik(4) . Agreement  is  to  within  a few 
parts  in  103  over  a + 25°  angular  range,  fhe  curve 
appears  to  be  slightTy  shifted  about  the  e = 0 line. 

This  is  attributed  to  a small  misorientation  in  the 
crystal.  Departure  from  theory  over  tne  remainder  of 
the  angular  range  could  be  accounted  for  by  a slight 
misorientation  in  the  plane  of  propagation.  The 
object  distribution  chosen  for  this  experiment  was  a 
random  array  of  small  gold  dots  and  L was  approximately 
100  wavelengths.  Although  we  reached  the  expected 
accuracy  estimated  from  simple  th'-  retical  considerations, 
the  gold  dots  produced  a very  non-^-'iform  object 
spectrum.  We  believe  that  the  ratner  "noisy"  slowness 
curve  is  accounted  for  by  this  fact. 

Figure  3 shows  the  measured  slowness  curve  for  YZ 
lithium  Niobate  again  compared  with  the  theoretical 
curve  computed  by  Slobodmk(4) . The  object  distribution 
ir,  this  case  was  a 2\  transducer  which  enabled  us  to 
measure  the  slowness  curve  over  + 18°.  As  in  the 
previous  example,  a slight  shift  jn  the  slowness  curve 
about  the  0 = 0 line  is  observed.  L was  400  wavelengths 
and  the  estimated  accuracy  is  in  this  case  within  a few 
Wrts  in  104.  Although  the  noise  in  the  slowness  curve 
ha  been  considerably  reduced,  a discrepancy,  up  to 
G-3  at  180  is  observed.  Subsequent  sensitivity 
ralc>. ’ations  on  the  theoretically  computed  curve 
indicated  that  this  discrepancy  is  well  within  the 
1’mnts  of  uncertainties  arising  from  tne  limited  known 
accuracy  of  the  elastic  tensor  elements  and  the  orient- 
ation of  the  crystal . 


We  have  demonstrated  a technique  for  measuring 
the  slowness  curves  of  acoustic  surface  with  an 
accuracy  of  a few  parts  in  ID-*  - substantially  better 
than  can  be  attained  by  other  means.  The  results  are 
indeed  more  precise  than  the  known  values  of  the 
velocities  as  derived  from  separate  measurements  on 
elastic  and  piezoelectric  tensor  elements.  This 
suggests  the  possibility  of  using  o measured  slowness 
curve  and  "working  backwards"  to  improve  the  accuracy 
of  the  measured  elastic  constants.  The  basic  approach 
is  firstly  to  obtain  the  differential  effect  of  th- 
variation  of  each  independent  element  on  the  velocity 
curve,  followed  by  an  optimization  procedure  to  obtain 
the  closest  possible  fit  to  the  experimental  data.  For 
any  particular  cut,  the  effect  of  certain  elements 
will  predominate.  One  can  increase  both  accuracy  and 
the  number  of  elements  measured  ty  undertaking 
iieasurements  on  several  cuts.  In  a very  similar  way 
we  can  also  derive  the  crystal  orientation.  It  seems 
probable  that  this  technique  will  permit  an  accuracy 
s.’bstc  ntially  better  than  the  + J°  normally  obtained 
f.'jm  X-ray  mease rements. 
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In  the  experiments  described  so  far,  the  trans- 
fers were  fired  along  a symmetry  axis  of  the  crystal. 

. ''Seqj»r;tly,  the  measured  slowness  cu'-ves  were  all 
'metrical  about  the  •<  = 0 line.  Clearly  one  could  make 
•*tter  use  of  the  available  object  spectrum  b\  firing 
transducers  at  an  angle  to  tne  symmetry  axis,  tnere- 
\>  deriving  slowness  data  over  a wider  angular  range. 
"'i*re  4 shows  tne  measured  slowness  curve  for  YZ 
'’tnium  Niobate  obtained  by  firing  tne  2<  transducer 
*’  angle  of  200  to  the  Z axis.  All  other  parameters 
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Figure  1:  Configuration  for  slowness  surface 
measurement 


Figure  3:  Velocity  results  for  YZ  Lithium  Niobate 
compared  with  Slobodnik  et  al. 
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Figure  2:  Velocity  results  for  YX  Quartz  compared 
with  Slobodnik  et  al. 


Figure  4:  Velocity  results  for  Lithium  Niobate 

obtained  using  a narrow  aperture  trans- 
ducer firing  at  20°  to  the  z axis. 
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BACKGROUND  SUBTRACTION  IN 
SURFACE-WAVE  HOLOGRAPHY 

tfuicuny  terms  At  • u'ln  \urtaie  Reau-onh  \'»irajr*r 

Suri.tu*-*j\c  holographic  technique*  can  he  used  io  ddcu  the 
presence  of  sm.iij  tie-fee; . 'I  lie  detection  limit  is  usual!)  deter- 
mined not  h>  signal  to  noise  considerations,  hut  rjthcr  b>  the 
presence  ol  spunous  background  vgruK.  such  as  hulk  waxes. 
I he  letter  describes  a technique  for  eliminating  this  back- 
ground. and  demonstrates  the  improvement  in  the  detection 
limit. 

Defect'!  on  a surface  can  be  identified  by  observing  their 
influence  on  the  propagation  of  acoustic  surface  waves.  A 
particular  technique  for  implementing  this  procedure  is  the 
use  of  2-dimensional  acoustic  holography.1, 1 With  reference 
u big.  I.  the  surface  is  illuminated  b>  acoustic  surface  waves. 
t-s-Mg  an  mteruigital  transducer  for  a pie/iielectric  material, 
or.  m other  cases,  by  a wedge  transducer.  The  complex  field 
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Fig.  1 Complex  held  distribution  is  measured  along  line  A A' 

distribution  n measured  along  the  line  AA\  this  record 
contains,  in  principle,  complete  information  on  the  field 
distribution  over  the  whole  surface  up  to  the  transmitting 
transducer,  provided  only  that  the  defects  arc  small  and 
sparse.  One  can  regard  the  record  as  a I -dimensional  holo- 
gram from  which  we  can  reconstruct  any  line  on  the 
2-dimensional  surface. 


The  distribution  can  be  recorded  on  film  and  subsequently 
reconstructed  optically.1  More  conveniently,  one  can 
digitise  the  record  and  reconstruct  using  purely  numerical 
techniques.'  It  is  ‘mportant  to  appreciate  that  the  latge-scale 
computations  needed  for  the  reconstruction  of  typicai 
2-dimensional  holograms  become  quite  insignificant  for  a 
I -dimensional  record.  Using  256  samples  and  8 bit  accuracy, 
the  reconstruction  can  be  effected  in  seconds  on  a large  com- 
puter and  in  times  of  the  order  of  a minute  on  a typical 
laboratory  minicomputer. 

The  area  of  search  for  defects  will  normally  lie  in  the  near 
field  of  the  transmitting  transducer.  The  field  distribution, 
even  in  the  absence  of  defects,  is  therefore  complex.  More- 
over, in  practice,  the  record  will  be  degraded  by  reflections 
from  the  edge  of  the  sample,  and  by  the  presence  of  bulk 
waves.  In  consequence,  there  is  a distinct  limit  to  the  smallest 
defect  which  can  be  discerned,  which  is  independent  of  the 
signal/noisc  ratio  attained  in  the  recording.  This  limitation 
can  be  overcome  if  one  has  available  a defect-free  sample  for 
comparison.  The  process  of  reconstruction  is.  of  course, 
linear;  one  is  therefore  at  liberty  to  subtract  the  -perfect’ 
record  from  that  obtained  on  the  test  sample  and  reconstruct 
the  difference,  h is  important  to  appreciate  that  the  reference 
sample  need  not  in  fact  lie  totally  free  of  defects  All  that  is 
needed  is  that  there  shall  be  so  few  that  one  can  discount  the 
possibility  of  a defect  on  the  reference  sample  in  the  same 
position  as  one  on  the  test  sample. 

To  demonstrate  this  technique,  we  adopted  the  configura- 
tion of  Fig.  1 . using  a /-cut  quartz  sample.  The  complex  field 
distribution  was  measured  by  a laser  probe  system, J the 
separate  amplitude  and  phase  records  being  digitised  and 
stored.  The  reconst  action  process  consists  essentially  of 
expanding  the  field  >.ito  a plane-wave  spectrum,  multiplying 
each  elementary  wave  by  the  appropriate  phase  delay,  allow- 
ing for  the  anisotropy  of  the  material,  and  subsequently  re- 
combining at  the  desired  reconstruction  i»;..ne. 

The  experiment  was  conducted  at  60  MHz  with  a ‘defect’ 
consisting  of  a gold  dot  approximately  30  //m  in  diameter  and 
4000A  thick,  t.e.  uppioximatcly  A/2  in  diameter  and  /.,'IOO  m 
thickness.  The  ‘perfect’  comparison  sample  was  obtained 
very  simp!)  by  .‘iping  off  the  gold  dot.  Fig.  2 a shows  the 
hologram  with  the  defect  present  and  Fig.  2b  shows  its 
reconstruction  m the  plane  of  the  dot.  It  would  clcarlv  be 
impossible  to  deduce  the  presence  of  a defect  from  thiv 
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evidence  The  hologram  with  the  gold  dol  removed.  Fig.  2c, 
in  not  easily  distinguished.  by  eye.  from  that  of  I ig.  la.  How- 
ever. the  reconstruction  of  the  difference  between  them. 
Fig.  hi.  shows  the  presence  of  the  defect  clearly.  The  back- 
ground icvel  in  Fig.  2 </  is  due  to  long-term  instabilities  in  the 
laser  probe  system  rather  than  U>  the  basic  noise  level.  Further 
improvements  in  the  detection  limit  can  therefore  be  expected. 

The  use  oi  the  .nine  sample  for  both  the  reference  and  the 
test  object  simplified  the  demonstration.  In  practice,  one 
would  have  to  take  additions  steps  to  ensure  that  the  location 
ot  the  recording  line  corresponds  as  closely  as  inissiblc  in  the 
two  cxpcniiients.  A lateral  displacement,  i.c.  along  the 
recording  line,  can  Ik*  icadiiy  eliminated  hy  progressively 
shifting  one  tetord  with  respect  to  the  other,  to  obtain 
maximum  correlation.  Similarly,  small  displacements  of  the 
recording  line  in  the  v -direction  can  Ik*  allowed  for  by  record- 
ing three  records  it*  span  the  uncertainty  in  the  .v-position. 
and  again  correlating  to  obtain  the  best  (it.  The  amount  of 
computation  needed  for  they  adjustments  is  again  insigni- 
ficant 

The  preliminary  experiments  encourage  the  view  that  the 
method  m.n  be  capable  of  development  into  an  instrument 
capable  of  detecting  extremely  small  defects  on  a surface  by 
tecording  one,  or  possibly  three,  linear  surface-wave  field 
disir’hiilions. 
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small  - though  artificial  - defects.  There  is  still  a good  deal  of  room  for 
further  refinement  of  the  method.  However  the  work  has  been  confined  so  far  to 
piezo-electric  materials.  It  is  clear  that  the  next  stage  must  be  to  extend  the 
work  to  non-piezo  electric  materials.  Amongst  those  that  are  most  likely  to  be 
of  direct  interest  are  glasses  and  semiconductor  materials.  The  main  problem 
is  to  devise  an  effective  surface  wave  transducter.  It  was  hoped  that  results  on 
this  - though  obtained  after  the  termination  of  the  contract  - would  still  become 
available  in  time  to  include  in  the  final  report.  This  unfortunately  was  not  possible 
- as  a result  of  a number  of  purely  technological  problems.  We  are  however 
continuing  to  work  towards  this  end. 

The  work  on  slowness  determination  is,  in  our  view,  also  of  possible 

relevance  to  n.d.t.  Of  course  slowness  is  a measurement  averaged  over  a surface. 

It  is  not  therefore  able  to  identify  individual  defects.  However  it  is  possible 

to  identify  quite  small  changes  in  the  average  characteristics  of  a surface  - for 

example  as  a result  of  surface  stress.  There  are  also  possibilities  of  measuring 

occlusion  of  gases  in  semiconductor  slices  after  processing.  The  current  precision 

4 

in  velocity  which  we  have  achieved  is  about  one  or  two  parts  in  10  . We  are  however 
confident  that  with  further  effort  we  should  be  able  to  achieve  at  least  one  addi- 
tional order  of  magnitude  in  precision.  We  believe  that  in  this  technique 
we  are  on  the  verge  of  developing  an  extremely  sensitive  surface  measurement  device. 
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C OVERALL  ASSESSMENT  - PLANS  FOR  FURTHER  WORK 

The  major  portion  of  the  effort*  on  this  contract  has  been  deployed  on  the 
microwave  scanning  technique.  The  results  which  we  have  obtained  are,  we  believe, 
sufficient  to  encourage  the  further  pursuit  of  this  method.  We  have  been  able 
to  identify  and  locate  cracks  in  metals  with  dimensions  of  one  micron  or  less. 

We  believe  therefore  that  there  may  be  direct  applications  for  surface  defect 
detection  in  metals.  Beyond  that  we  have  demonstrated  the  use  of  the  method 
•for  inspection  of  electrical  parameters  of  materials  - notably  dielectric  constants, 
permeability  ferrimagnetic  resonance  effects  etc.  In  this  area  we  believe  that 
the  instrument  we  have  developed  is  capable  of  performing  not  only  a competitive 
but  a unique  task.  We  know  of  no  other  way  in  which  the  surface  value  of. 
these  parameters  can  be  established  with  comparable  accuracy. 

There  are  a number  of  ways  in  which  the  instrument  could  be  developed  so  as 
to  lead  to  a more  fully  engineered  version,  - and  giving  superior  performance. 
Briefly  the  improvements  one  would  seek  include  the  following: 

(i).  Use  of  higher  frequency.  Using  available  solid  state  sources  one 
might  work  at  8mm  wavelengths. 

(ii)  Use  of  a pair  of  integrated  resonators,  one  for  measurement  the 
other  for  carrier  cancellation.  This  would  improve  stability. 

(iii)  Simplification  of  electronic  system;  use  of  lower  noise  components 
which  have  recently  become  available. 

(iv)  Improved  mechanical  scanning  stage. 

The  work  we  have  done,  has  been  exploratory.  We  have  tried  to  assess  the 
attainable  performance  for  a wide  range  of  quite  different  materials  testing 
problems.  In  our  view,  further  development  should  probably  be  concentrated  on  a 
more  limited  set  of  objectives,  and  aimed  at  attaining  the  highest  performance 
within  that  set. 


*The  doctoral  student  whose  salary  was  paid  on  the  contract  worked  on  chi s 
exclusively. 
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The  acoustic  holography  and  velocity  determination  work  has  also  reached 
a stage  where  we  believe  serious  consideration  should  be  given  to  developing 
the  technique  for  n.d.t.  purposes.  The  main  bottleneck  has  been  the  lack  of 
effective  transducers  for  non-piezo  electric  materials  - which  of  course  represent 
the  largest  class  of  samples  one  would  consider  testing.  Here  again,  we  feel 
that  further  work  should  probably  be  concentrated  on  a limited  class  of  materials 
to  be  tested.  We  believe  that  glasses  and  semiconductor  materials  are  likely  to 
prove  particularly  suitable.  In  the  case  of  glasses  one  would  be  able  to  detect 

microcracks  and  also  stress.  In  the  case  of  semiconductors  the  sensitivity  of 

the  technique  has  now  reached  a point  where  quite  small  defects  will  show  up. 

Again  the  overall  velocity  measurement  is  believed  now  to  be  sufficiently  sensitive 
to  allow  the  assessment  or  a number  of  surface  characteristics  - including 
stress,  and  distribution  of  impurities.  There  is  a possibility  that  one  might 
be  able  to  assess  the  effect  of  occluded  gases  as  well. 

We  are  intending  to  submit  a new  proposal  for  the  extension  of  the  acoustic 
surface  wave  work. 


